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quadratic control problem lies n . scheme in the paper approximates 

Riccati operator equation. The appro * ^ sequence of finite dimensional LQG 

the infinite dimensional LQG problemh ‘dimensional , usually finite element 
problems defined for a mo del of the structure. Two 
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The finite dimensional e q uations estimator gains 

veloped , including formulas for co mparison of gains based on dit- 

to their functional representation to alluw^ of the approximating control and 

ferent orders of approximation. d ine finite dimensional compensators is 

estimator gains and of the corresponding finite systeras produced 

studied. Also, convergence and sta Y discussed. The convergence 

with the finite dimensional compensat - f he finlte dimensional 

theory is based on the convergence of the solutions Riccatl equa _ 

tions^ A "nume r i c a ^ °e x'anq) 1 e° w i t h° a flexible beam, a rotating rigid body, and a 
lumped mass is given. 
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1. introduction 


Recent yeans have seen increasing research in active control of flexible 
structures. The primary motivation for this research is contr 
flexible aerospace structures, which are becoming larger and more flexible at 
th e same time that their performance regents are becoming more stringent, 
for example, in tracRing and other applications, satellites with large enten- 
te solar collectors and other flexible comments must perform fast slew 
maneuvers while maintaining tight control over the vibrations of their flexi- 
ble elements. Both of these inflicting objectives can be achieved only with 
a sophisticated controUer. There are applications also to control of robotic 
manipulators with flexible lints, and possibly to stabilization of large civil 
engineering structures such as long bridges and tall buildings. 

The first question that must be answered when designing a controller for 
a flexible structure is whether a finite dimensional model is sufficient as a 
hasis for a controller that will produce the required performance, or is a 
distributed model necessary! While some structures can he modeled well by a 
fixed n«t*r of dominant modes, there are structures whose flexible character 
nan * captured sufficiently for precise control only by a distributed model. 
Still others - perhaps most of the aerospace structures of the future - 
he modeled sufficiently for control purposes by some finite dimensional 
approximation, but an edentate approximation may he impossible to determine 
before design of the stroller, or .mPensator. This paper deals with struc- 
tures that are flexible enough to require a distributed model in the design o 


an optimal LQG compensator. 


The linear-quadratie-gaussian optimal control problm, for distributed, or 
xnfiMte dimensional. system is a Generalisation to Hilbert space or the now 
olaaaical WO problem for finite dimensional systems. Tbe solution to the 
infinite dimensional problm, yields an infinite dimensional stat^estimator- 
based compensator, which is optimal in the context of this paper. By a 
separation principle IB,. 041, tt e problem reduces to a detepmlnistic linear 
Quadratic optimal control probl- and an optimal estimation, or filtering 
problem with gauaaian white noise. In infinite dimensions, the control system 
dynamics are represented by a semigroup of bounded linear operators instead of 
«.e matrix exponential orators in finite dimensions, and the plant noise 
process may be an infinite dimensional random poooess. solutions to both 

the control and filtering problems Involve Hloeati operator motions, which 
-e generalisations of the Hiccati matrix equations in the finite dimensional 
Current results on the infinite dimensional LOG problem are most com- 
plete for problems where the input and measurement operators are bounded, as 
this paper requires throughout. This boundedness also peemits the strongest 

approximation results here, for related control problmss with unbounded input 
and measurement, see [C3, C5 , LI, L2] . 


our primary objective in this pa*r is to approximate the optimal infi„ 
ite dimensional U0 compensator for a districted model of a flexible struc- 
ture with finite dimensional compensators Csed on approximations to the 
structure, and to have these finite dimensional compensators produce near 
optimal performance of the closed-loop system. We discuss how the gains that 
the finite dimensional compensators converge to the gains that 
determine the infinite dimensional compensator, and we examine the sense in 
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which the finite dimensional compensators converge to the infinite dimensional 
compensator. With this analysis, we can predict the performance of the 
closed-loop system consisting of the distributed plant and a finite dimen- 
sional compensator that approximates the infinite dimensional compensator. 

Our design philosophy is to let the convergence of the finite dimensional 
compensators indicate the order of the compensator that is required to produce 
the desired performance of the structure. The two main factors that govern 
rate of convergence are the desired performance (e.g. , fast response) and the 
structural damping. We should note that any one of our compensators whose 
order is not sufficient to approximate the infinite dimensional compensator 
closely will not in general be the optimal compensator of that fixed order; 
i.e. , the optimal fixed— order compensator that would be constructed with the 
design philosophy in [B7 , B8] . But as we increase the order of approximation 
to obtain convergence, our finite dimensional compensators become essentially 
identical to the compensator that is optimal over compensators of all orders. 

An important question, of course, is how large a finite dimensional com- 
pensator we must use to approximate the infinite dimensional compensator. In 
[G6, G7 , G8, Ml], we have found that our complete design strategy yields com- 
pensators of reasonable size for distributed models of complex space struc- 
tures. This strategy in general requires two steps to obtain an implementable 
compensator that is essentially identical to the optimal infinite dimensional 
compensator; the first step determines the optimal compensator by letting the 
finite dimensional compensators converge to it; the second step reduces, if 
possible, the order of a large (converged) approximation to the optimal com- 
pensator. The first step, which is the one involving control theory and 
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approximation theory for distributed systems, is the subject of this paper. 

For the second step, a simple modal truncation of the large compensator some- 
times is sufficient, but there are more sophisticated methods in finite dimen- 
sional control theory for order reduction. For example [G8, Ml], we have 
found balanced realizations [M2] to work well for reducing large compensators. 

The approximation theory in this paper follows from the application of 
approximation results in [B 6, G3 , G4] to a sequence of finite dimensional 
optimal LQG problems based on a Ritz-Galerkin approximation of the flexible 
structure. For the optimal linear^quadratic control problem, the approxima- 
tion theory here is a substantial improvement over that in [Gl] because here 
we allow rigid-body modes, more general structural damping (including damping 
in the boundary), and much more general finite element approximations. These 
generalizations are necessary to accommodate common features of complex space 
structures and the most useful finite element schemes. For example, we write 
the equations for constructing the approximating control and estimator gains 
and finite dimensional compensators in terms of matrices that are built 
directly from typical mass, stiffness and damping matrices for flexible struc- 
tures, along with actuator influence matrices and measurement matrices. 

For the estimator problem, this paper presents the first rigorous approx- 
imation theory. (We have used less complete versions of the results in previ- 
ous research [G 6, G7 , G8, Ml].) As in the finite dimensional case, the infin- 
ite dimensional optimal estimation problem is the dual of the infinite dimen- 
sional optimal control problem, and the solutions to both problems have the 
same structure. Because we exploit this duality to obtain the approximation 
theory for the estimation problem from the approximation theory for the 
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optimal control problem, the analysis In this paper is almost entirely deter- 
ministic. We discuss the stochastic interpretation of the estimation problem 
and the approximating state estimators briefly, but we are concerned mainly 
with deterministic questions about the structure and convergence of approxima- 
tions to an infinite dimensional compensator and the performance - especially 
atability - of the closed-loop systems produced by the approximating compel 

sators. 

next, an outline of the paper should nelp. The paper has two main parts, 
which correspond roughly to the separation of the optimal UO problem into an 
optimal linear-quadratic regulator problem and an optimal state estimation 
problem. The first half. Sections 2 through 6. deal with the control system 
sod the optimal regulator problem. Sections 7 throw*, 10 treat the state 
estimator and the compensator that is formed by applying the control law of 
the first half of the paper to the output of the estimator. 

Section 2 defines the abstract model of a flexible structure and the 
energy spaces to be used throughout the paper. We assume a finite number of 
actuators, since this is the case in all applications, and we assmne that the 
actuator influence operator is bounded. Section 2 also establishes oertai 
mathematical properties of the operloop system that are useful in control and 
approximation. To our pledge, the exponential stability theorem in Section 
2.3 is a new result, and we find it interesting that such a simple Lyapunov 
functional accommodates such a general damping model. 

Section 3 discusses the linear-quadratic optimal control problem for the 
distributed model of the structure and establishes some estimates involving 
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bound, on solutions to Infinite dimensional Rlccatl equations and open-loop 
and closed-loop decay rates. We need these estimates for the subsequent 
approximation theory. To get the approximation theory for the estimation 
problem, we have to give certain results on the control problem in a more gen- 
eral form than would be necessary were we Interested only in the control prob- 
lem for flexible structures. Therefore, in Section 3. as in Sections 5 and 7. 
»e first give some generic results applicable to the WG problem for a variety 

of distributed systems and then apply the generic results to the control of 
flexible structures. 

Because we assume a finite number of actuators and a bounded input opera- 
tor. the optimal feedback control law consists of a finite number of bounded 
linear functionals on the state space, which is a Hilbert space. This means 
that that the feedtank law can be represented in terms of a finite number of 
vectors, which we call functional control gains, whose inner products with the 
generalized displacement and velocity vectors define the control law. For any 
finite- rank, bounded linear feedback law for a control system on a Hilbert 
space, the existence of such gains is obvious and well known, a functional 
control gain for a flexible structure will have one or more distributed com- 
ponents. or kernels, corresponding to each distributed component of the struc- 
ture and scalar components corresponding to each rigid component of the struc- 
ture. 


We introduce the functional control gains at the end of Section 3. and we 
introduce analogous functional estimator gains in Section 7. The functional 
sains play a prominent roll in our analysis. They give a concrete representa- 
tion of the infinite dimensional compensator and provide a criterion for con- 
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vergence of the approximating finite dimensional compensators. 

We develop the approximation scheme for the control problem in Section 4. 
The idea is to solve a finite dimensional linear-quadratic regulator problem 
for each of a sequence of Ritz-Galerkin approximations to the structure. We 
develop the approximation of the structure in Section 4.1 and prove conver- 
gence of the approximating open-loop systems, The approximation scheme 
includes most finite element approximations of flexible structures. For con- 
vergence, we use the Trotter-Kato semigroup approximation theorem, which was 
used in optimal open-loop control of hereditary systems in [B5] and has been 
used in optimal feedback control of hereditary, hyperbolic and parabolic sys- 
tems in [B6, G1 , G3] and other papers. The usual way to invoke Trotter-Kato 
is to prove that the resolvents of the approximating semigroup generators con- 
verge strongly. To prove this, we introduce an inner product that involves 
both the strain-energy inner product and the damping functional, and show that 
the resolvent of each finite dimensional semigroup generator is the projec- 
tion, with respect to this special inner product, of the resolvent of the ori- 
ginal semigroup generator onto the approximation subspace. The idea works as 
well for the adjoints of the resolvents, and when the open- loop semigroup gen- 
erator has compact resolvent, it follows from our projection that the approxi- 
mating resolvent operators converge in norm. 

In Section 4.2, we define the sequence of finite dimensional optimal con- 
trol problems, whose solutions approximate the solution to the infinite dimen- 
sional problem of Section 3. The solution to each finite dimensional problem 
is based on the solution to a Riccati matrix equation, and we give formulas 
for using the solution to the Riccati matrix equation to compute approxima- 
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tions to the functional control gains as linear combinations of the basis vec 
tors. 

Section 5.1 summarizes some generic convergence results from [B6, G3 , G4] 
on approximation of solutions to infinite dimensional Riccati equations. Sec- 
tion 5.2 applies these generic results to obtain sufficient conditions for 
convergence — and nonconvergence — of the solutions of the approximating 
optimal control problems in Section 4.2. A main sufficient condition for con- 
vergence is that the structure have damping, however small, that makes all 
elastic vibrations of the open-loop system exponentially stable. This is a 
necessary condition if the state weighting operator in the control problem is 

coercive. 

In Section 6, we present an example in which the structure consists of an 
Euler- Bernoulli beam attached on one end to a rotating rigid hub and on the 
other end to a lumped mass. We emphasize the fact that we do not solve, or 
even write down, the coupled partial and ordinary differential equations of 
motion. For both the definition and numerical solution of the problem, only 
the kinetic and strain energy functionals and a dissipation functional for the 
damping are required. We show the approximating functional control gains 
obtained by using a standard finite element approximation of the beam, and we 
discuss the effect on convergence of structural damping and of the ratio of 
state weighting to control weighting in the performance index. As suggested 
by a theorem in Section 5 , the functional gains do not converge when no struc- 
tural damping is modeled. 

In Section 7, we begin the theory for closing the loop on the control 
system. We assume a finite number of bounded linear measurements and 
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construct the optimal state estimator, which is infinite dimensional in gen- 
gj. 32 . The gains for this estimator are obtained from the solution to an 
infinite dimensional Riccati equation that has the same form as the infinite 
dimensional Riccati equation in the control problem. We call these gains 
functional estimator gains because they are vectors in the state space. 

Since the approximation issues that this paper treats are fundamentally 
deterministic, we make the paper self contained by defining the infinite 
dimensional estimator as an observer, although the only justification for cal- 
ling this estimator and the corresponding compensator optimal is their 
interpretation in the context of stochastic estimation and control. We dis- 
cuss the stochastic interpretation but do not use it. We say estimator and 
observer interchangeably to emphasize the deterministic definition of the 
estimator here. 

With the optimal control law of Section 3 and the optimal estimator of 
Section 7, we construct the optimal compensator, which also is infinite dimen- 
sional in general. The transfer function of this compensator is irrational, 
but it is still an m(number of actuators) p(number of sensors) matrix func- 
tion of a complex variable, as in finite dimensional control theory. The 
optimal closed-loop system consists of the distributed model of the structure 
controlled by the optimal compensator. 

Approximation of the optimal compensator is based on approximating the 
infinite dimensional estimator with the sequence of finite dimensional estima- 
tors defined in Section 8.1. The gains for the approximating estimators are 
given in terms of the solutions to finite dimensional Riccati equations that 
approximate the infinite dimensional Riccati equation in Section 7. Although 
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defined as observers, these finite dimensional estimators can be interpreted 
as Kalman filters, as shown in Section 8.2. In Section 8.3, we give formulas 
for finite dimensional functional estimator gains that approximate the func- 
tional estimator gains of Section 7. These approximating estimator gains 
indicate how closely the finite dimensional estimators approximate the infin- 
ite dimensional estimator. In Section 8.4, we apply the Riccati equation 
approximation theory of Section 5 to describe the convergence of the finite 
dimensional estimators. 


Most of the results in Section 8 are analogous to results for the control 
problem and follow from the same basic approximation theory, but certain 
differences require careful analysis. There is an important difference in the 
way that the Riccati matrices to be computed are defined in terms of the fin- 
ite dimensional Riccati operators. Indeed, the Riccati matrix equations to be 
solved numerically might seem incorrect at first. To demonstrate that the 
finite dimensional estimators that we define in Section 8.1 are natural 
approximations to the optimal infinite dimensional estimator, we show in Sec- 
tion 8.2 that each finite dimensional estimator is a Kalman filter for the 
corresponding finite element approximation of the flexible structure. The 
brief discussion in Section 8.2 is the only place in the paper where stochas- 
tic estimation theory is necessary, and none of the analysis in the rest of 
the paper depends on this discussion. 

In Section 9.1, we apply the nth control law of Section 4 to the output 
of the nth estimator to form the nth compensator. (The order of approximation 
is n.) The nth closed-loop system consists of the distributed model of the 
structure controlled by the nth compensator. Since each finite dimensional 
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estimator is realizable, the nth compensator and the nth closed loop system 
are realizable. In Section 9.2, we discuss how the sequence of realizable 
closed-loop systems approximates the optimal closed-loop system. Probably the 
most important question here is whether exponential stability of the optimal 
closed-loop system implies exponential stability of the nth closed-loop system 
for n sufficiently large. We have been able to prove this only when the 
approximation basis vectors are the natural- modes of undamped free vibration 
and these modes are not coupled by structural damping. That this stability 
result can be generalized is suggested by the results in Section 9.3, which 
describe how the transfer functions of the finite dimensional compensators 
approximate the transfer function of the optimal compensator. 

In Section 10, we complete the compensator design for the example in Sec- 
tion 6. Assuming that white noise corrupts the single measurement and that 
distributed white noise disturbs the structure, we compute the gains for the 
finite dimensional estimators and show the functional estimator gains. As in 
the control problan, the functional gains do not converge when no damping is 
modeled. We apply the control laws computed in Section 6 to the output of the 
estimators in Section 10 to construct the finite dimensional compensators, and 
we show the frequency response of these compensators. As predicted by Section 
9.3, the frequency response of the nth compensator converges to the frequency 
response of the optimal infinite dimensional compensator as n increases. In 
Section 10.3, we discuss the structure and dimension of the finite dimensional 
compensator that should be implemented. 

We conclude in Section 11 by discussing where the approximation theory 
presented in this paper is most complete and what further results would be 
most important. 
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2. The Control System 


We consider the system 

x (t) + D 0 x(t) + A 0 x(t) = B 0 u(t). t > 0, 

( 2 . 1 ) 

where x(t) is in a real Hilbert space H and u(t) is in R m for some finite m. 
The linear stiffness operator A Q is densely defined and selfadjoint with com- 
pact resolvent and at most a finite number of negative eigenvalues. We will 
postpone discussion of the damping operator D Q momentarily, except to say that 

it is symmetric and nonnegative. The input operator B Q is a linear operator 
from R 01 to H, hence bounded. 

By na tural 32 ties , we Jiill sear} tile eigenvectors £f XM eigenvalnA 
problem 


Vj = Vj • 

( 2 . 2 ) 

From our hypotheses on A Q , W e know that these eigenvalues form an infinitely 
increasing sequence of real numbers, of which all but a finite number are 
positive. Also, the corresponding eigenvectors are complete in H and satisfy 

<4> V ^j > H = <A 0^i' *j>H = °» 1 * J- 

(2.3) 

(These properties of the eigenvalue problem (2.2) are standard. See, for 

example, [Bl] , [Kl] .) For A, . > 0 id - -VT i a „ 

J r is a natural frequency . 
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Remark 2._1. Our analysis includes the system 
• • • 

^Ux(t) + D n x(t) + A n x(t) = B ft u(t) , t > 0, 

U U U ( 2.1 ) 

where the mass operator Mq is a selfadjoint, bounded and coercive linear 

operator on a real Hilbert space H Q . The operators A Q , B Q and D Q in (2.1') 

have the same properties with respect to Hq that the corresponding operators 

in (2.1) have with respect to H. To include (2.1') in our analysis, we need 

only take H to be Hq with the norm- equivalent inner product <*, = 

, and multiply (2.1') on the left by In H, the operator J^Aq 

is selfadjoint with compact resolvent, and Mq^Dq is symmetric and nonnegative. 
With no loss of generality then, we will refer henceforth only to (2.1) and 
assume that the H-inner product accounts for the the mass distribution. 0 

2.1 The Energy Spaces and the First-order Form of the System 

The Elastic-Strain - Energy Space V and Total - Energy Space E 


We choose a bounded, selfadjoint linear operator A^ on H such that Aq = 
Aq + A^ is coercive; i.e. , there exists p > 0 for which 


<AqX, x>h 2. P 1 1 x | | H » 


x e I (A 0 ) = D(A 0 ). 


(2.4) 


Since Aq is bounded from below, there will be infinitely many such A^'s. 
In applications like our example in Section 6, it is natural to select for A^ 
an operator whose null space is the orthogonal complement (in H) of the eigen- 
space of Aq corresponding to nonpositive eigenvalues. Obviously, any A^ that 
makes Aq coercive must be positive definite on the nonpositive eigenspace of 
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With Aj chosen, we define the Hilbert space V to be the completion of 
D(A q ) with respect to the inner product < v i» v 2^v = <AqV^ » v 2 >h» and 

<vi / 2 / 2 ' s 'i / 2 

V 2 c D(Aq) . Note that V = D(Aq ) and <v^,V 2 >y = <Aq v^.A^ v 2^h* (Since A^ 
is a bounded operator on H, different choices of A^ yield V’s with equivalent 
norms, thus containing the same elements). 

In the usual way, we will use the imbedding 


VC H = H' C V' , 


where the injections from V into H and from H into V* are continuous with 
dense ranges. We denote by A y the Riesz map from V onto its dual V'; i.e. , 


< v,v l>V = <Vi>v, v l' v e V * 

Then A^ is the restriction of Ay to D(Aq) in the sense that 

(Ay v i)v = <v,AqV 1 > h , v x e D(A 0 ), v e V. 


(2.5) 


( 2 . 6 ) 


Now we define the total energy space E = V x H, noting that when A^ is 

• 2 

coercive and x(t) is the solution to (2.1), then | | (x(t) ,x(t))| |- is twice the 
total energy (kinetic plus potential) in the system. We want to write (2.1) 
as a first-order evolution equation on E. To to this, we must determine the 
appropriate semigroup generator for the open-loop system. We will derive this 
generator by constructing its inverse explicitly, and then we will try to con- 
vince the reader that we do have the appropriate open-loop semigroup genera- 
tor. The approach seems mathematically efficient, and we will need the 
inverse of the generator for the approximation scheme. First, we must state 
our precise hypotheses on damping and discuss its representation. 
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Actually, we do require an operator D Q defined from some subset of H 
into H. Rather, we assume only that there exists a damping functional 

d 0 (v l' v 2 ): v * V-*R 

(2.7) 

such that d Q is bilinear, symmetric, continuous (on V X V) and nonnegative. 

If we have a symmetric, nonnegative damping operator D Q defined on D(A 0 ) 
such that D Q i s bounded relative to A 0 , then <D 0 v 1 .v 2 > H defines a bilinear, 
symmetric, bounded, nonnegative functional on a dense subset of V X V. In 
this case, the unique extension of this functional to V X V is d Q . (That D Q 

being A Q -bounded implies continuity of <Dq‘,*> h with respect to the V-norm 
follows from [Kl, Theorem 4.12, p. 292].) 


Under our hypotheses on dQ, there is a unique bounded linear operator Aq 
from V into V such that 


V v ’V = (A d v 1 )v, v 1 ,v e V. 

( 2 . 8 ) 

The operator (Ay Aq) is then a bounded linear operator from V to V, and 
My Aq) is self ad j oint (on V) because dQ is symmetric. Also 

d 0 (v.v 1 ) = <v,Ay 1 A D v 1 >y = <Ay 1 A D v.v 1 >y v,,v e V. 

' (2.9) 

R emark 2.2. We chose to begin our description of the control system model 
with (2.1) because its form is familiar in the context of flexible structures. 
The stiffness operator Aq, for example, is the infinite dimensional analogue 
of the stiffness matrix in finite dimensional structural analysis. In appli- 
cations like the example in Section 6, though it is often easier to begin 
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with a strain-energy functional from which the correct strain-energy inner 
product for V is obvious. The stiffness operator is defined then in terms of 
the Riesz map for V (see [S3] for this approach), rather than V being defined 
in terms of the stiffness operator; specifically, X Q is defined by (2.6) with 
D(Xq) = Ay^H. Either way, the relationship between and V is the same. But 
the only thing that needs to be computed in applications is the V-inner pro- 
duct; an explicit Aq need not be written down. O 


The Semigroup Generator 


We define A * e L(E,E) by 


r-i 




( 2 . 10 ) 


This operator is clearly one-to-one, and its range is dense, since V is dense 


in H and D(Aq) is dense in V. Now, we take 


x = a-h- 1 . 


( 2 . 11 ) 


Direct calculation of the inner product shows 


<A h ' h ^E d Q (h,h) , 


( 2 . 12 ) 


so that X is dissipative with dense domain. Also, since D(A ) = E, A is max 
imal dissipative by [Gl, Theorem 2.1]. Therefore, X generates a C Q - 
contraction semigroup on E. 


Finally, the open-loop semigroup generator is 
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/W 

A = A + 


0 ol 

A 1 °J 


D(A) = D(A) , 


(2.13) 


where A^ is the bounded linear operator discussed above. With 


B = 


r° ] 

KJ 


e L(R m ,E). 


(2.14) 


the first-order form of (2.1) is 


;(t) = Az(t) + Bu(t), t>0 


(2.15) 


where z = (x,x) e E. 


To 3 ee that A is indeed the appropriate open-loop semigroup generator, 
suppose that Aq is coercive (so that A^ = 0) and that we have a symmetric, 
nonnegative Aq— bounded damping operator Dq. Then the appropriate generator 
should be a maximal dissipative extension of the operator 


o 

A 



D(A q ) X D(Aq) . 


(2.16) 


o 

It is shown in [Gl, Section 2] that A has a unique maximal dissipative exten- 
sion, and it can be shown easily that the A defined above is an extension of 
after noting that, in the present case. 


(A v Ap) | 


D(A 0 ) 



(2.17) 


We should note that Showalter [S3, Chapter VI] elegantly derives a 
semigroup generator for a class of second— order systems that includes the 
flexible-structure model here. The presentation here is most useful for our 
approximation theory because of the explicit construction of the inverse of 
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the semigroup generator. For the purposes of this paper, we do not need to 
characterize the operator A itself more explicitly, but we should make the 
following points. 


First, from A we see 


D(A ) = { (x, x) : x e V, x + A^ApX c D(A Q ) } 


(2.18) 


In many applications, especially those involving beams, the "natural boundary 
conditions" can be determined from (2.18) and the boundary conditions included 
in the definition of DAq) . in the case of a damping operator that is bounded 
relative to , D(A) = D(Aq)xV. If the damping operator is bounded rela- 
tive to Ag for n < 1, then A has compact resolvent. 


In mar^y structural applications, the open-loop semigroup is analytic, 
although this has been proved only for certain important cases. Showalter 
obtains an analytic semigroup when the damping functional is V-coercive; for 
example, when there exists a damping operator D Q that is both Abounded and 
as strong as A Q . Such a damping operator results from the Kelvin- Voigt 
viscoelastic material model. Also, it can be shown that the semigroup is ana- 
lytic for a damping operator equal to c Q Ag for 1/2 < (i 1 1 and o fl a positive 
scalar. The case n = 1/2, which produces the same damping ratio in all modes, 
is especially common in structural models, and Chen and Russell [Cl] have 
shown that the semigroup is analytic for a more general class of damping 
operators involving A*^. 


Finally, we can guarantee that the open— loop semigroup generator is a 
spectral operator (i.e. , its eigenvectors are complete in E) only for a damp- 
ing operator that is a linear combination of an H-bounded operator and a 
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fractional power of A^. However, nowhere do we use or assume anything about 
the eigenvectors of either the open-loop or the closed-loop semigroup genera- 
tor. The natural modes — of undamped free vibration — in (2.2) are always 
complete in both H and V. 

2.2 The Adjoint of 1 

i y — ♦ 

Since (Ay A^) is selfadjoint on V, direct calculation shows that A = 
(A -1 )* — the adjoint of *A _1 with respect to the E- inner product — is 

~ A V A D ^ 1 

= l -! 0 J * 

(2.19) 

Then 1* = Ck~* ) -1 . Having K - * explicitly facilitates proving strong conver- 
gence for approximating adjoint semigroups* 

2.3 Exponential Stability 

The following theorem says that, if there are no rigid-body modes and if 
the damping is coercive (basically, all structural components have positive 
damping), then the open-loop system is uniformly exponentially stable. That 
the decay rate given depends only on the lower bound for the stiffness opera- 
tor and the upper and lower bounds for the damping functional is essential for 
convergence results for the approximating optimal control problems of subse- 
quent sections. The theorem is a generalization of Theorem 6.1 of [Gl] to 
allow more general damping, but the proof is entirely different and much 
nicer. The current proof uses an explicit Lyapunov functional for the homo- 
geneous part of the system in (2.15). Recall that T(’) is the open-loop semi- 
group, with generator A, and E is the total energy space V x H. 
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Theorem 2.3 . Suppose that A^ and dQ are H-coercive. Let p be the positive 
constant in (2.4) , and let 6^ and 8^ be positive constants such that 

8 oII v Hh 1 d o (v * v) ^ 8 ili v ilr v 6 v * 

( 2 . * 

Then 


I |T(t ) 1 | 1 (1 + 5 1 & 0 /2) 1/2 exp[-t/(^- +A+ Sj) J , t 2 0. 

° ( 2 . 21 ) 

Proof . For y > max{ ^ , jp ) , define Q e L(E) as 


Q = 


(yl+Ay Ap) Aq 

i yi 


( 2 . 22 ) 


Since (Ay*Ap) is selfadjoint and nonnegative on V, Q is selfadjoint and coer- 
cive on E. Define the functional p(*) on E by 


p(z) = <Qz, z> E = r(||x||y + ||x||g) + 2<x,x> h + d Q (x,x), 

(2.23) 

* 

where z = (x,x). From (2.4), we have 

2|<x,x> h | < | lx | | V J |x| | H 1 yjyl |z | |g, (2 24) 

so that 

<r-^)||z||* i p(s) i t-'uzul (2 _ 2J) 


with 


Next take a - (*.i> . D<A, and set A(x,i» - <*»> • E ’ ° r 

(x.i) - A _1 <y.y) - <-Vv-*o ly - y) ' 

Note x = y e V. Now, 

<qaz , z> E - <Q(y-y>' A_1(y ' y)> E 
- *<y.A;V>v - HAv'W y| ‘v + ,|yl, H 
= - tllxllj + yV X,X) " ,,X,I H 3 * 


(2.27) 


(2.28) 


From (2.20) then. 


<QAz,z> E <. 


- tllxllv + (rV 1 ),,x,, H 1 L 3 3z ' 


Therefore p (*> is a Lyapunov 
(2.26) and (2.29) , with y 


functional, and the theorem 



. O 


follows from 


(2.29) 
(2.25) , 


3. The Optimal Control Probl 


em 


ection 3.1 presents some preliminary def initio 

optimal linear-quadratic regulator problem, on an arbitrary 1 T ^ 

These results are generic in e-h Albert space. 

the sense that the Hilbert space El. f 

— - — *. - - 0 „ ratora , B . et : - 

~ — structure as ln 

r hair of ^ ^ ~ ~ — 13 , m aiw us ; o :: 

— - the infinite dimensional atate J 

analogous resui ts for the control probl«, Subs „ ^ ^ 016 

prooiem. Subsection 3.2 gives * 

tant implications of the general results for t„ 

^ th6 CQSG whpno f l A 

is that defined In Section 2 . control system 


3-1 tte Generic Optimal Segulator Problem 


1st a linear operator 4 generate a C 

-os K. and suppose B . Q . M ^ ^ ™ — 

w e M£ " E ) and R s L(R m ) n 

- aelfadjoint and » positive deflate and symmetric ihe ml " 

*. E is to choose the control u . t (0 
functional 2 ( '° >iR ^ to “inimi^e the cost 


J(z(0) , u) 


/ (<Qz(t) ,z(t)> E + 


<nu(t).u(t)> R m)dt. 


where the state 2 (t) is given by 


( 3 . 1 ) 


z(t) 


T(t)z(0) + 


t 

/T(t-rj) Bu(r) )d R t t 20 . 


£efinltJ,bn i.i. 


( 3 . 2 ) 


A function u e L 2 ( 0 .-; D ) 


iS an .control f or ^he 
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initial state z, or simply an admissible control for z, if J(z,u) is finite; 
i.e. , if the state z(t) corresponding to the control u(t) and the initial con- 
dition z(0) = z is in L 2 (0,«°;E). 

Definition 3 .2 Let the operators A, B, Q, and R be as defined above. An 
operator II in L(E) is a solution of the Riccati algebraic equation if II maps 
the domain of A into the domain of A* and satisfies the Riccati algebraic 

equation 

A*n + ha - n br - 1 b* n + q = o. 

(3 .3) 

Theorem 3.3 (Theorems 4.6 and 4.11 of IG43 ) . There exists a nonnegative sel- 
fadjoint solution of the Riccati algebraic equation if and only if, for each 
z e E, there is an admissible control for the initial state z. If II is the 
minimal nonnegative selfadjoint solution of (3.3), then the unique control 
u( ) which minimizes J(z,u) and the corresponding optimal trajectory z( ) are 
given by 

u ( t ) = -R _ 1 B* nz(t) 

(3 .4) 

and 

z(t) - Sit)*, (3 , 5) 

where S(') is the semigroup generated by A-BR 1 B^I. Also, 

J(z,u) = min J(z,v) = <IIz,z>g 

v (3.6) 

If, for each initial state and admissible control. 
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11m | |z(t) 1 1 =0, 

t— >°° 

(3.7 

there exists at most one nonnegative selfadjoint solution of (3.3). If Q is 
coercive, (3.7) holds for each initial state and admissible control and S( • ) 
is uniformly exponentially stable. □ 


We will refer to T(-) as the open - loop semigroup and to S(») as the 
QPt Iwal closed-loop semigroup . 


To prepare for the convergence analysis in Sections 5 and 9, we must 
present now some rather arcane estimates for the decay rate of the closed-loop 
system in the optimal control problem. 

T ^gpreq 3. .4 Suppose that the open-loop semigroup T(») satisfies 
I lT(t ) 1 1 1 M 1 e° lt , t 1 0, 

(3.8) 

for positive constants ^ and oj, that H is the minimal nonnegative selfad- 
joint solution to (3.3), and that S(t) is the optimal closed— loop semigroup in 
Theorem 3.3. If there exists a constant such that, for each z e E, 

OD 

f IIS(t)z|| 2 dt 1 ^(<nz,z> E + 1 1 z 1 1 2 ) 

(3.9) 

and a constant ^ such that 

II n ii < 

(3.10) 

then there exist positive constants ft, and a 2 , which are functions of , 1^, 
M^, and only, such that 
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Ils(t)|| < v^e 1 , t 1 o. 


(3.11) 


Proof . This follows easily from Theorems 2.2 and 4.7 of [ ]. □ 

kenfflg 3 .5 . Suppose that there exist positive constants M and a such that 

I |T(t) 1 1 1 Me' at , t 1 0. 

t 

If z(0) 8 E, h e L o (0,®;E), and z(t) = T(t)z(0) + / T(t-s)h(s)ds, then 
1 0 


/ | | z ( t ) | | 2 dt 1 
0 

Proof . The result follows from (3.12), the convolution theorem [Dl, page 95l] 
and the triangle inequality. □ 

tomma 3.6. Suppose that E is finite dimensional and that the pair (Q,A) is 
observable (in the usual finite dimensional sense). Then there exists a con- 
stant M, which is a function of A, B and Q only, such that 

/| |z(t) | | 2 dt L M (/ (<Qz(t) ,z(t)> E + | | u ( t ) 1 1 2 ) dt , 

0 0 (3.14) 

where z(t) is given by (3.2). 

Proof . The proof, which is at most a mild challenge, is based on the fact 
that the observability grammian 
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(3.15) 


W(t)= / e ATt Qe At dt 

0 

is ooercive for any positive t. □ 


The next theorem says, among other things, that if the open - loop control 
system decouples into a finite dimensional part that is stabilizable (in the 
usual finite dimensional sense) and an infinite dimensional part that is uni- 
formly exponentially stable, then the entire system is uniformly exponentially 
stabilizable, so that (3.3) has a nonnegative selfadjoint solution. 


X toprep) 3. .7,. Suppose that there exists a finite dimensional subspace EqCT 
D(A) such that E Q and reduce A (and T(t)), and write 


where A^ 
Simil arly. 


r A n o ] 

r B n] „ 

■ Q ll 

Q 12' 

1 0 A 22j * 

l b 21 J- q ■ 

. °12 

<$2. 


and A 22 are the restrictions of A to E Q and D(A) X Eg, 


(3.16) 

respectively. 


T(t ) 


-Tn<t) 

0 


° i. 


(3.17) 


Also, suppose that the pair ( A 11 ,B 11 ) is stabilizable and that there exist 
positive constants M^, and p such that 


I I ^22 ( t ) | | 


-a^t 


t 2 0 


and 


(3.18) 
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max{| | B | |,| | Q| | } i p. 


(3.19) 


i) Then there exists F e L (E, H 111 ) such that A-BF generates a uniformly 
exponentially stable semigroup on E. Also, (3.3) has a nonnegative selfad- 

t 

joint solution, and the minimal such solution satisfies (3.10) with Mq a func- 
tion of Aj^, Bj^, R, M^, and P only. 

ii) If Q ^2 “ 0 and the observable, then there exists a 

unique nonnegative selfadjoint solution II to (3.3), and there exist positive 

9 9 

constants and Oj — which depend on A^, B^, Oj^* °i an< * ^ only 

such that the optimal closed-loop semigroup satisfies 

l|S(t)|| < Mje" 2 , t 2 0. (3 . 20) 

Proof , i) To say that (A^,b^) is stabilizable means that there exists a 
linear operator F n from E Q to R m such that each eigenvalue of A 11 -B n F n has 
negative real part. Hence A-BF generates a uniformly exponentially stable 
semigroup if F = [F^j 0] , so that there exists an admissible control for each 
initial condition in E. 

It is easy to write down an upper bound for the performance index in 
(3.1) in terms of R, P , M^, <* x and the decay rate of exp( [A n -B 1X F X1 It). 
That the in (3.10) depends only on A^, B 11# R, ^ and p follows then 
from the fact that F^ is a function of A^j and B^. 

9 9 

ii) Clearly, (3.8) holds with and depending only on k^, B^, 

and F^. Therefore, we have (3.8) and (3.10) with the bounds depending only 

on A n , B n , R, M^, a' and p . Finally, the existence of an % for 
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(3.9) which depends only on these parameters follows from using (3.1) and 
(3.6) in applying Lemma 3.6 to the part of the system on E Q and then Lemma 3.5 

to the part of the system on E^. Part ii) of this theorem then follows from 
Theorem 3 .4 . □ 

itefflarK 3.8 When we say in Theorem 3.4 that Mj and a 2 are functions of Mq» i^, 
Mj^ and oj only, we mean, for example, that for two optimal control problems 
on different spaces E, with different operators A, B, etc., if the same con- 

9 

stants and o x work in (3.8)-(3.10) for both problems, then the 

same constants M 2 and a 2 will work in (3.11) for both problems. Similarly, in 
Theorem 3.7 ii) , as long as Ep, A^, Bj^, R» and j) remain the 

same, the same Mj and a 2 will work in (3.20) even if Aj 2 , B n and Q 22 
change. Q 

3.2 Application to Optimal Control of Flexible Structures 

For the rest of this section, Ap, Aj, A. T(t). Bp and B are the operators 
defined in Section 2,1, and E = V X H is the energy space defined there. 

Remark 3.9. Theorem 3.7 is useful mainly when all but a finite number of 
modes have coercive damping in the open- loop system and the damped and 
undamped parts of the open-loop system remain orthogonal. This is the case, 
for example, with modal damping. The next theorem does not require ortho- 
gonality of the damped and undamped parts of the system, but it does require 
an independent actuator for each undamped mode. The situation of Theorem 3 .10 
is typical in aerospace structures: Any elastic component should have some 

structural damping, but rigid-body modes are common; for a structure to be 
controllable, an actuator is required for each rigid-body mode. □ 
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Theorem 3 .10. i) 

Aj are H-coercive, 
that, for all v e 

1 1 v | |y 2 Pllvllj. 
cf Q (v,v) 2 Pllvllj. 
d Q (v.v) 1 y||v||y* 

and 

max{||B 0 ||.||Q|M|R||} i P 

(The V-continuity of d Q implies (3.23).) 
selfadjoint solution II, which satisfies 
and P only . 

ii) Suppose also that 


r** 

A^^ and d Q = d Q + 

, y and p such 

( 3 . 21 ) 
( 3 . 22 ) 

(3.23) 

(3.24) 

Then (3.3) has a minimal nonnegative 
(3.10) with a function of p, t 

(3.25) 

(3.26) 

p , Y and p only . 


<0z,z> E 2 pi |z| lg* z e E. 

Then the optimal closed-loop semigroup satisfies 
||S(t)|| i V “ 2 . t 2 0. 

where Mj and o 2 are positive constants depending on 
Proof , i) The suboptimal control 


* 

Suppose that A^ = BqBq and that A^ = Aq + 
so that there exist positive constants p 
V, 


2 l ) 



u(t) - -B*[x(t) + x ( t ) ] 

(3.27) 

produces a closed-loop system with exponential decay at least as fast as that 
in Theorem 2.1. The required upper bound in (3.10) follows then from (3.1), 
(3.6) and (3.24) 

xi) In this case, the ^ in (3.9) is p , and Theorem 3 .4 yields the result> g 


Now we will consider the structure of the optimal control law in 


more 


n 


where 


e L(E 

,E) 

Do 

% 

* 

n i 

“2 

, n x e 

! L( 1 

With 

z = 


(3.28) 


u ( t ) = -R n*x(t) + I^xU)]. 


(3.29) 


Since B q e L(R m ,H), we must have vectors b ± * H, li i 1 m , such that 


m 


B n u = Z b, u 


i=l 


ii 


(3.30) 


for 


U " ( V 2 *" 


(3.31) 


Also, for h e H, 


B 0 h - [<b 1 ,h> H <b 2 ,h> R <b m ,h> H ] T . 

Since HjXtt) and I^xtt) are elements of H, we see from (3.29) and 


(3.32) 
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(3.32) that the components of the optimal control have the feedback form 


where f^ 


We call 


u^t) — (f^.x(t))y x(t i 1. • • . »m, 


(3.33) 


e V and g^ e H are given by 


1 < R_1 > ij n i b J, 


(3.34a) 


m 

I 

j=l 


(R 1 ) 


ij 


% b J. 


i = 1, 


t m. 


^ and g^ functional gains . 


(3.34b) 
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4. The Approximation Scheme 


4.1 Approximation of the Open- loop System 


Hypothesis 4.J,. There exists a sequence of finite dimensional subspaces of 
V such that the sequence of orthogonal projections Py R converges V-strongly to 
the identity, where P yn is the V-projection onto V n . Also, each V n is the 
span of n linearly independent vectors e.. 

J 

Since it should cause no confusion, we will omit the subscript n and write 
just ej, keeping in mind that the basis vectors may change from one V R to 
another, as in most finite element schemes. Also, we will refer to the Hil- 
bert space E n = X V R , which has the same inner product as E = V X H. 

For n 2 1* we approximate x(t) by 


n 

x(t) = I ?,( t)e,, 
j=l J J 

where £(t) = (£ x (t), £ 2 (t), ... 5 n (t)) T satisfies 

M n 5(t) + D n $(t) + K^t) = Bju(t), 


(4.1) 


(4.2) 


and the mass matrix M n , damping matrix D n , stiffness matrix K n , and actuator 
influence matrix Bq are given by 
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D n = * 


M n = [<e i ,6j> H ] , 


k" - uJ'V^Vh 1 " [<e i’ e J ! V> - KVi'Vh 1, 


(4.3) 


B 0 “ t <e i* b j > H 1 * 


Of course, (4.2) can be written as 
; = k\ * B n u 


( 4 . 2 # ) 


where 


and 


n - u T . i 1 i T 


(4.4) 


A n = r ° 1 1 , B n = 

[-M" n K n -M n D n J 


M-Xl 


(4.5) 


Note; 


Throughout this paper, we use the superscript n in the designa 
matrices In the .*•> approximating system and control prohl^, U* 

A”, b", H n , etc. Hence the superscript n Indicates the order of ^ 
approximation - and it never indicates a power of the matrix. By H ”. 
we denote the inverse of the mass matrix M 

in the designation of a linear operator in the n‘» approximation, -e use the 
the subscript n. For example. A n and B„ are the operators whose matrix 
representations are A n and B n . respectively. 
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For convergence analysis. It is useful to note that (4.1, and (4.2, or 
(4.2') are equivalent to 


V fc) = Vn (t) + B n u(t), 


wh6r*6 z =: ( y y \ n , . 

n <VV V and 4 n , L(E n , and B n , L(B ". are the ^ 
whose matrix representations are given in ,4.5,. Also, for any real A. 


(4.6) 



(4.7) 


and 


uV 


+ AD“ ♦ it ",. 1 - <AM n + D n ,p 2 + m"(I 2 


(4.8) 


if 


o 2 = Xa 1 - 0 1 


(4.9) 


V J - r j i 

n ~ ^ aie and h J 

i=l 1 1 n 


J = 1 * 2 . 


(Substituting A" and (4.10, into (4.7, yields <4.«) and (4.9)). 


(4.10) 


Next, we will prepare to invoke the Trotter-Kai-n 

ccer-Kato semigroup approximation 

theorem to show how (4.2), (4 2') ann (a <\ 

' ’ ) and (4.6) approximate (2.1) and (2.15) 

First, ue ,111 treat the case in ubioh a„ 1s coercive ,„o rigid-body modes,, 
so that A,-0 and A^A,,; the general case Is a straight-forward extension. For 
h„ coercive, the -p.-l.op semigroup generator A is maximal dissipative. Also, 

fnM An nk _ A 


For each n, A„ is dissipative on F, main idea here is to projeot (A-A^ 


-1 
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onto V n in a certain inner product and observe that the result is exactly 
(X-A n ) _1 , where A r is the operator on V n in (4.6) and (4.7). Of course, we 

need only do this for real X>0. 

For real X>0 then, define an innerprcduct on V by 


2 y . • < 


<*,*>, = X < * » * > H + ^0 


xd n c,‘) + <*»*> 


V* 


x H U ’ ( 4 . 

Under the hypotheses in Section 2 on d 0 , is norm-equivalent to 

For n 2 1. let P n (x) be the projection of V onto V n in the inner product 
<*,*> Now let h 1 ,h 2 e H and note that 


(>-*) [ V 2 


(4 


is equivalent to 

(:*)■ *'1$) ■($ • 

With A -1 from (2.10), (4.13) is equivalent to 

(I ♦ xA« l Ap ♦ xV "' 1 - + NF’V * 1 * A » lh2 (4 


and 


v 2 = Xv 1 - h 1 . 


(' 


If 


v 1 = P (X) v 1 
n n 


and v 


n 


P n (X)v‘. 


it follows from (4.11) and (4.14) that 


11 ) 

V* 


. 12 ) 

.13) 

\ .14) 

) .15) 
4.16) 


3:» 





= * ^i'W^H + ^iVV^V + < e i'V 1 >; 

= <e i ,a 2 A 0 " 1 + W- 1 A D+ I)v 1 > v 

= <e i (XAo 1 My 1 A D )h 1 +A- 1 h 2 > v , 
and from (4.15) that 


< e i'V >* 


(4.17) 


<e i' v l\ = <e i* v2> x = >-<e i ,v 1 > 


X <e i’ h \ . 


(4.18) 


»«. for h* . h l . v n . h 2 . h 2 „ . v n . and vj. vj, hj and h 2 written aa in 
(4.10), (4.17) and (4.18) yield (4.8) and (4.9) again. 

This shows that 


r p n ( x) i 

l 0 P n (X)J a_A) 1 |fi n = a_A n ) -1 * 


(4.19) 


which yields 


r p n a> 

0 ] 

l » 

P n (A) J 


^En - ^ _A -) Xp 


n' En 


(4.20) 


-hare P En is 11 s IrfinUfifiiiin f into »e projection P En can be written 

P . P Vn 0 | 

En " L » p HnJ • 


(4.21) 


where P„ n is the V-projeotlon onto ea before, and P„ n ia the H-projectlon 
onto V n . Since the V-nom la stronger than the H-norn and the norm induced, by 
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the X-innerproduct is equivalent to the V-norm, it follows from Hypothesis 4.1 
that U-A n ) 1 p converges E-strongly to (X-A) -1 as n-» «. Now, with A 

I 11 

extended to E n as, say. n(P En -I) , Trotter-Kato [EL, page 504, Theorem 2.16] 
yields the following. 

Theorem 4 . .2,. For Aq coercive, let T^ (*) be the (contraction) semigroup gen- 
erated on E r by A R . Then, for each t 1 0, T n (t)Pg n converges strongly to 
T(t), uniformly in t for t in bounded intervals. 

In the general case, when A Q i s no t coercive, the open-loop generator A 
is obtained from the dissipative £ by the bounded perturbation in (2.13), so 
that [G3 , Theorem 6.6] yields the following generalization of Theorem 4.2. 

gorollary 4.3. Let T n C) be the semigroup generated on E r by A n . Then, for 
each t 2 0, T Q (t )P En converges strongly to T(t), uniformly in t for t in 
bounded intervals* 

4 - 4 - When A has compact resolvent, (X-A n ) -1 P En converges in L(E) to 

(X-A) -1 . 

iCSifif. This follows from (4.20) and a standard result that the projections of 
a compact linear operator onto a sequence of subspaces converge in norm if the 
projections converge strongly to the identity, a 3 do P En and P n (X) . □ 

That the adjoint semigroups also converge strongly follows from an argu- 
ment entirely analogous to the proof of Theorem 4.2. In particular, equations 
like (4.11)-(4.17) are used to show that 
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(4.22) 


[ P n<*> 0 1 

l 0 P n<*>J 


(X-A*) _1 P 


En 


" (x A n ) P En 


In showing this. A * is used as A -1 was used above. Also, care must be taken 
to calculate A* with respect to the E-inner product. The result is 


Theorem 4.5. Let T (*) be the sequence of semigroups in Corollary 4.3. Then, 
for each t 1 0, T*(t)P En converges strongly to T (t), uniformly in t for t in 
bounded intervals. 


Finally, for the approximation to the actuator influence operator B e 
L(R m » E) , recall B n e L(R m > E n ) , the operator whose matrix representation is 
the matrix B n in (4.5). From (4.3), it follows that 

B n = P En B * (4.23) 

Since B has finite rank m, B n and B* converge in norm to B and B*. respec- 
tively. 
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4.2 The Approximating Optimal Control Problems 

The n th acmaJ control pr obl SB 13 : « iven V 0) ‘ (x n (0) '*n< 0) > e E n’ 

choose u e L 2 <0 fi" 1 ) to minimize 


J n (z„(°). u) 


"(<Q n z„<t>- z n (t)> E + <»u(t).u(t)> B m)dt, ^ ^ 


where Q n - P En QI En - *« a33ume: 

. \ i t mu p there exists an admissible con 

Hypothesis 4.6. For each n > 1 and z n (0)e E n . there exi 

trol (Definition 3.1) for (4.6) and (4.24). □ 

.... Pnr> Hypothesis 4.6 is that, for each n, the system 
A sufficient condition for Hypothesis 

(A ,B ) be stabilizable. 

n’ n' 


By Theorem 


3.1, the optimal control u R (t) has the feedback form 


u n U> * - B " lB n n n z n <t) 


(4.25) 


. „ F n is nonnegative and selfadjoint, and 1^ 
where 1^ is a linear operator on E n , !!„ is nonne fe 

satisfies the Riccati equation 


«> n * Hn*n - * Q n ' 


(4.26) 

AS a result of Hypothesis 4.6. (4.26) has at least one nonnegative, selfad- 
joint solution. The minimal such solution is the correct ^ here. If the 
system (A„.Q n ) is observable, then n„ ia the unique nonnegative, selfadjoint 
solution to (4.26) , and is tositive definite. If »e write II,, as 
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\ - 


P^On Uln 


*ln 


then (4.25) becomes 


(4.27) 


Vt) = - r" B 0 [nJ n x n (t) + i^ n ; n(t)] . 

The feedback law (4.28) can be written in functional 
as in Section 3. We have 


(4.28) 

-feedback form, just 


where 


u n (t) - t u in<t> u 2n (t) --V t,]T ' 


(4.29) 


and 


U in (t) <f in' x n (t)> V “ ^in'V^H’ 1 < i i m . 


(4.30) 


m 


in 


g in 


^ ^R ^i n PHn^j' T i. i <. m. 


j=l 


m 


^2n P Hn b j ' 1 1 i 1 m. 


course, f^ and are the n th approximations to the functional 
and g A in (3.34) . 


(4.31a) 


(4.31b) 


gains f 


^ction 5, (4.25) (4.31) will be useful for studying how the solution 
to the n optimal control problem converges to the solution to the original 
problem of Section 3. but for numerical solution of the n th problem, »e need 
the matrix representations of these equations. 

We will need the following grammian matrices: 
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and 


(4.32) 


w n 



(4.33) 


Note: The matrix W n will be the inverse of W n . The superscript n on any matrix 

indicates the order of approximation, not a power of the matrix. Also, recall 
the note following (4.5) . 


Now recall = PE n Q|g n . Since Q = Q* e L(E) and E = V X H, we can write 


K 

where Q Q = oje L(V), Q 1 8 L(H,V), and Q 2 = Q 2 « L(H) . 
lation shows that 


(4.34) 

Straightforward calcu- 


Q n = w _n Q n ' 


(4.35) 


where Q n is the matrix representation of Q and Q n is the nonnegative, sym- 

fi. 

metric matrix 





with 


Q 


n 

0 


^ e i » ^0 e j ^ ' 


(4.36) 
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Q^ej)y] , 


(4.37) 


°2 [<e i'Q2 e j>H ] * 

Also, recall that A n and B n are the operators whose matrix representations are 
given by (4.5), and note that the matrix representations of A* and B* are 
W -n (A n ) T W n and (B n )^W n , respectively. 

With the matrix representation of denoted by if* the Riccati operator 
equation (4.26) is equivalent to the Riccati matrix equation 


W" n (A n ) T W n II n + Ifk n -D n B n R- 1 (B n ) T W n D n + Q n = 0. 


(4.38) 


While °n is self adjoint, if in general is not symmetric, but the matrix 


I f = W n lf 


(4.39) 


is symmetric and nonnegative, and positive definite if II n is. 
(4.39) by W n , we obtain 


Premul ti plying 


(A n ) T if 1 + I^A 11 - fi n B n R -1 ( B n ) T If 1 + Q n = 0, 


(4.40) 


which is the Riccati matrix equation to be solved numerically. 


Now we need one more set of matrix equations for the numerical solution 
of the n th optimal control problem. Since the functional gains f^ and g in 
are elements of V they can be written as 


n f. 


£ “j ®J “ d *1- 


n g ± 

1 p; , i = 1, . . . ,m, 

j=l J J l A A 


f i g± f i Si 

where p , p e R n . We need equations for p and p in terms of Hr. One 
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way to get these equations is to partition ^(obtained from (4.3 9)) and then 
work out the matrix representation of (4.31). However, another approach is 
more instructive because it relates the present Hilbert space methods to the 
standard finite dimensional solution of the n th optimal control problem. 

th 

The n optimal control problem can be stated equivalently as: given 

n(0) = U(0) T ,5(0) T ] e R 2n choose u e L 2 (o,<*;R m ) to minimize 


J n (ri(0),u) = J (il(t) T Q n Ti(t) + u(t) T Ru(t))dt, 

0 

where n(t) = U(t) T ,£ (t) T ] T satisfies (4.2'). For (4.2') and (4.42), 
optimal control law is 


(4.42) 

the 


% = -R-VnVt) 

(4.43) 

where IT 1 is the minimal nonnegative, symmetric solution to (4.40). 

Since S is related to x r by (4.1), the optimal control u n in (4.43) must 

be equal to the optimal control u n in (4.29) -(4.31) . Substituting (4.41) into 
(4,30) yields 


u in = - O - 0 8i )M = -KpS 1 O gi ) T ]W^. 

Then, using (4.44) and equating (4.29) and (4.43) yields 


(4.44) 


f l f 2 

P P 

e 2 

P P 2 *'- 


m 


a m 


= w^IPbHr -1 


(4.45) 


— -QSH Mve ihe . complete .aplqilon to the n th optimal control problem : 
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The Riccati matrix equation (4.40) is solved for DP; then the optimal control 
is given by (4.43). and equivalently by (4.29M4.30) with the functional 
gains f in and g in given by (4.41) and (4.45). In the next section, we will 
give sufficient conditions for the solution to the n th optimal control problem 
to converge to the solution to the optimal control problem in Section 3 for 
the original infinite dimensional system. 
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5. Convergence 


As in Section 3, subsection 5.1 will state some results for the optimal 
linear regulator problem involving generic linear operators A, B. Q. etc. , on 
an arbitrary real Hilbert space E, and subsection 5.2 will expand upon these 
results for the particular class of control problems treated in this paper. 

5.1 Generic Approximation Results 

Let the Hilbert space E and the linear operators A. T(*). B. Q and R be 
as in Section 3. Suppose that there is a sequence of finite dimensional sub- 
spaces E n . with the projection of E onto E„ denoted by P En . such that P En con- 
verges strongly to the identity asn4«, and suppose that there exist 

sequences of operators A r e L(E n > , B n e L(R m »E n >f Q n = Q n e L(E n ) ' Q n ^ 
such that we have the following strong convergence. For all z e E and t2P, 

exp(A n t)P En z -> T(t)z (5.1) 

and 

exp(A*t)P En z -4 T*(t)z (5.2) 

as n -4- «, uniformly in t for t in bounded intervals; for each u e R m , 

B n u -4 Bu; (5.3) 

for each z e E, 


Theore.. 5 .1 ■ Suppose that for each n there Is a nonnegati.e, aelfadjolnt 
linear operator!^ on E„ which eatiaflea the Rlooati aleebralo equation 
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0 . 


(5.5) 


+ V„ - W lB X * “ n * 

If there exist positive constants M and 0, independent of n, such that 
l|exp([A n -B n R _1 B*D ii ]t)|| < Me“ Pt , t l 0 , 

(5.6) 

and if 1 1 13^1 1 is bounded uniformly in n, then the Riccati algebraic equation 
(3.3) has a nonnegative selfadjoint solution II, and, for each z e E, 

Dn p En z Hz 

(5.7) 

and 

expUin-B^BXmp^z -► S(t)z 

(5.8) 

uniformly in t 2 0, where S(*) is the semigroup generated by A-BR _1 B^I . if 
there exists a positive constant 6, independent of n, such that 

Q n^ 8 ' 

(5.9) 

then || being bounded uniformly in n guarantees the existence of positive 

constants M and 0 for which (5.6) holds for all n. 

£rfiof. The theorem follows from Theorem 5.3 of [G4] when the operators A n , Q 

and Hj, are extended to all of E by defining them appropriately on e£. For the 

details of this procedure, see Section 4 of [Gl] . Or better. Banks and Kun- 

isch [B6] have modified Theorem 5.3 of [G4] to obtain essentially the present 

theorem without using the artificial, and rather clumsy, extensions to E‘ L i n 

n ±n 

the proof. □ 
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T he or en) 5.2 The strong convergence in (5 7 ) < mn i j. . „ 

K e in { 5 . 7 ) implies uniform norm convergence 

of the optimal feedback laws: 


^®n ^n^E “ 2 H|| o as n —> ®. 
n 


( 5 . 10 ) 


*“*■ ttiS f ° UC “ S fr °“ th9 SelfadJOl " t “- o': H and P En and the finite 
dimensionality of the control space See equations (4.23, and (4.24, of 
[Gl], □ 

asfirssS.l Assume the hypotheses of Iheorem 5.1 but do not assume (5.<S> or 
(5.9). If || ij^n la bounded uniformly in n. then the Ricoatl algebraic 
equation (3.3, has a nonnegative selfadjoint solution n. and. for each a e E 
Vito* converges weakly to Hz. 

inar. This is Theorem 4.7 of [031. whose proof is valid under the hypothese, 
here. Q 


Th ° aaonleoming of the weak convergence in Theorem, 5.3 
not yield uniform norm convergence of the feedback control 


is that it does 
laws. 
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5.2 Convergence of the Approximating Optimal Control Probleme of 
Section 4 .2 

For the rent of thin section, A,,, 1 1 , A. T(t). B 0 and B are the operators 
defined in Section 2. The operators A„. B„. Q n andl^ are the operators in 
the approximation schmne of Section 4. In particular. U . is the 

minimal nonnegative, self-adjoint solution of the Riccati operator eouation 
(4.24) . According to Corollas 4.3 and Theory 4.5. the Ritn-GalerRin approx- 
imation scheme presented in Section 4.1 converges as required in (5.1) and 
(5.2); (5.3) and (5.4) follow frcm (4.23) and the definition 0 n - P En Ql En in 
Section 4.2. Also. Hypothesis 4.4 guarantees for each n the existence of the 
required solution of the Riccati equation (5.5) in Theorem 5.1. 

Since n„ is nonnegative and self-adjoint, its eigenvalues, which are 
also the eigenvalues of its matrix representation, are real and nonnegative, 
and its norm is equal to its maximum eigenvalue. 

Theorem i.i If 0 is E-ooercive and d 0 - 0 (i.e.. there is no open-loop damp- 
ing), then there is no nonnegative self-adjoint solution of the Riccati opera- 

tor equation (3.3) , and 


|| njl • as n ». 


(5.11) 


proof . Recall the operator 'k in Section 2.1. By Theorem 1 of [02] . there can 

he no compact operator C a L(E.E) such that T e C generates a uniformly 
exponentially stable semigroup. Therefore, since a compact linear perturba- 
tion of t yields A. there can be no compact linear C such that A + C generates 
a uniformly exponentially stable semigroup. 


48 



Now, unless (5.11) holds, there exists a subsequence such that || ]T || 

J 

is bounded in nj, so that Theorem 5.3 says that there exists a nonnegative 
self-adjoint solution II of (3.3). Since Q is coercive. Theorem 3.3 then 
says that the semigroup generated by A - BR~ 3 B* II is uniformly exponentially 
stable. But this is impossible — BR -1 B* II is compact because its rank is not 
greater than m. □ 


Theorem 5 .5 . Suppose that Aq and are both H coercive. Then there 

exist positive constants and independent of n# such that 

-a, t 

||exp[A n t]|| 1 M^e , t 2 0. 


Proof . First, we define Aq r and Dq r e L(V n ,V Q ) to be the operators whose 
matrix representations are M _n K n and M _n D n , respectively. (See Section 4.1.) 
The operator A R is then 


A 


n 




(5.13) 


Since Aq and d Q are H-coercive, there exists a positive constant p, 
independent of n, such that 


<A 0n h ' h> H 


2 Pllhlln 


(5.14) 


and 

<D 0n h ' h> H ^ Pi I* 1 ! Ih (5 

for all h e V . Since d n is continuous on V X V, there exists a positive 
n u 
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constant y, independent of n, such that 


<D 0n h,h> H i yllhllj 

(5.16) 

for all h e V Q . The theorem follows then from Theorem 2.1. □ 

T heQr^m £.6 . Suppose that A Q has an invariant subspace V Q which is also 
invariant under the damping map A^Ap , that E Q = V Q X V Q is a stabilizable 
subspace for the control system, and that the restrictions of A Q and d 0 (»,») 
to Vq are both H-coercive. Also, suppose that V Q has finite dimension n Q and 
that, for each n 2 n Q in the approximation scheme, the first nQ e^ span V Q 
and the rest are orthogonal to V Q in both V and H. 

i) Then (3.3) has a nonnegative solution II, and for each n 2 Hq, (5.5) has a 
nonnegative self-adjoint solution 1^. Also, 1^ is bounded uniformly in n, so 
that converges toll weakly, as in Theorem 5.3. 

ii) If Eq and Eq (the E-orthogonal complement of Eq) are invariant under Q, 
and if the part of the open-loop system on E Q is observable with the measure- 
ment Qz, then (5.6)— (5.8) hold as in Theorem 5.1. 

Proof. We will invoke Theorem 3.7 to establish the existence of the uniform 
bounds and decay rates needed in Theorem 5.1. In the approximating optimal 
control problems, the part of the control system on E Q i s the same for each n; 
the approximation of the control system takes place on E^ We can write E = 

o n 

E 0 ® ^0n’ where E 0n is the ortho S° na l complement of E Q in E n> and E Q and E^ 
clearly reduce the open- loop semigroup for each n 2 Hq. 

For the part of the open- loop system on E Qn , Theorem 5.5 establishes 
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positive and a^, independent of n, for (3.17). Also, we have a (3 
independent of n for (3.18) because B n = P En B and Q n = P En Ql En * Therefore, i) 
follows from Theorem 3.7 i) and Theorem 5.3. 

The definition of Q n and the requirement on where the various basis vec- 
tors must lie imply that Eq and Ej^ n reduce Q n if Eq and E ^ reduce Q and that 
the restriction of Q to E Q is the same for all n. Therefore, ii) follows from 
Theorem 3.7 ii) and Theorem 5.1. □ 

Remark 5.7. In applications, the subspace V Q in Theorem 5.6 usually contains 
rigid-body modes. The theorem includes the case where both Aq and dQ are H- 
coercive on all of V (no rigid-body modes and all modes damped) . In this 
case, Vq is the trivial subspace. O 

Bamarif 5.8. Otherwise, for applications to flexible structures. Theorem 5.6 
usually requires two things: first, modal damping must be modeled for the 

structure, so that the natural modes remain uncoupled in the open-loop system; 
second, the natural mode shapes must be used for the basis functions in the 
approximating optimal control problems. Although these requirements may seem 
restrictive from a mathematical standpoint, such modeling and approximation 
predominate in engineering practice. Also, we get our strongest convergence 
results under these conditions. For applications where the basis vectors are 
not the natural mode shapes, the following theorem is useful. □ 

Theorem 5.9. Suppose that A Q + B q Bq and d Q + B q Bq are H-coercive. Then (3.3) 
has nonnegative solution II, for each n (5.5) has a nonnegative self-adjoint 
solution 11^, and || 1 1 is bounded uniformly in n. Hence Theorem 5.3 

applies. Furthermore, if Q is E— coercive, then (5.6)— (5.9) hold in Theorem 5 . 1 . 


51 



Proof . The required bounds follow from Theorem 3.10 and the proof of Theorem 
5.5. Although we took A^ = BqB^ in Theorem 3.10, this is not necessary in the 
final result, since all bounded self-adjoint operators A^ on H that make Aq + 
A^ coercive yield equivalent norms for V. 0 


Theorem 5.10. If (5.7) holds for each z e E, then 


Ilf 


in 


^lly °' 


(5.17a) 


H g in “ s i * Ih ~ ^ 0 ' as n 


(5.17b) 


where f^ and g^ are the functional gains in (3.16), and f^ n and g^ n are the 
approximating functional gains in (4.31) and (4.41) . 


Proof . The result follows from (4.31). □ 


Note that (5.10) and (5.17) are equivalent. 
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6 * Exam pi e 


6.1 The Control System 

One end of the uniform Euler- Bernoulli beam in Figure 6.1 is attached 
rigidly (cantilevered) to a rigid hub (disc) which is free to rotate about its 
center, point 0. which is fixed. Also, a point mass m, is attached to the 

other end of the beam. The control is a torque u applied to the disc, and all 
notion is in the plane. 



Figure 6.1, Control System 
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r = hub radius 
H = beam length 

Iq = hub moment of inertia about axis 
perpendicular to page through 0 

m^ = beam mass per unit length 
m^ = tip mass 

El = product of elastic modulus and second moment 
of cross section for beam 
fundamental frequency of undamped structure 


10 in 
100 in 

100 slug in 2 

.01 slug/in 
1 slug 

13 ,333 slg in 3 / sec 2 
.9672 rad/sec 


Table 6.1 Structural Data 


The angle 0 represents the rotation of the disc (the rigid-body mode), 
w(t,s) is the elastic deflection of the beam from the rigid-body position, and 
w^t) is the displacement of from the rigid-body position. For technical 
reasons, we do not yet impose the condition w^t) = w(t,|); more on this 

later. 

me control problem Is to stabilize rigid-body motions and linear (small) 
transverse elastic vibrations about the state e - 0 and v. - 0. Our linear 
model assumes not only that the elastic deflection of the beam is linear but 
also that the axial inertial force produced by the rigid-body angular velocity 
has negligible effect on the bending stiffness of the beam. The rigid-body 

angle need not be small* 

For this example, it 'is a straight forward exercise to derive the three 
coupled differential equations of motion in 0, w and , and they do have the 
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form (2.1'). However, to emphasize the fact that we do not use the explicit 
partial differential equations, we will not write these equations here. 

Rather, we will write only what is normally needed in applications: the 

kinetic and strain-energy functionals, the damping functional and the actuator 
influence operator. 


Remark 2.1 applies to this example, and to most examples with complex 
structures. The generalized displacement vector is 


x = (e»w,w 1 ) e H q = R L 2 (0,|)X R. 

The kinetic energy in the system is 

• 4 

Kinetic Energy = V 2 <x,x> H 

where H is Hq with the inner product 

a i 

<x,x> H = m b / Q [w+(r+s)6] tw+(r+s)©lds 


( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


+ + nijtwj+fr+Jhe] [Wj^+trn- )$] . 

As in most applications, we need not write the mass operator explicitly, but 
there exists a unique selfadjoint linear operator on Hq such that 

<x,x> H = <M 0 x,£> h . 

° (6.4) 

It is easy to see that Mq is bounded and coercive. Hence Hq and H have 
equivalent norms. 


The input operator for (2.1') (which maps R to H Q ) i s 
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( 1 , 0 , 0 ) . 


(6.5) 


Since we multiply (2.1') by to get (2.1), the input operator for (2.1) is 
(M^Bq) . Note that 

- Bj. < 6 . 6 

where (M^Bq)* 11 is the H-adjoint of (M Q 1 B 0 ) and B Q is the H Q -adjoint of B Q . 

Remark 2.2 also applies here. The only strain energy is in the beam and 
is given by 

Strain Energy = Vi a(x,x) ? 

with 

0 

a(x,x) = EI/ q w''w''ds, g 

where (*)" - 3 2 ( * ) /3s 2 ( * ) . To make aC,*) into an inner product, we must 
account for rigid- body rotation. Thus we set 

<x,x> v = a(x,x) +00 ^ 9 


and define 

V = tx = (0,0, 0(1) ) : 0 e H 2 (0,|), 0(0) = 0'(O) = 0 }. 
Also, we have 

<x,x> v = a(x,x) + <b o b o x '^ > H q 

= a(x,£) + < (Mq^Bq) ^0^0^ X,X ^H' 


( 6 . 10 ) 


( 6 . 11 ) 
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so that A x = B q B*, or (M^Bq) (Mq^) * H » depending on whether the H Q or the H- 
inner product is used in computing the V- inner product. But MS need neither 
nor A Q explicitly . We need only (£.8) and (6 .9) , along with (6 .3 ) , _fco compute 
the required inner products . 

As mentioned in Remark 2.2, the operator A Q can be defined now by (2.6), 
and the stiffness operator is A Q = A^ - Aj. Using the Hq- inner product in 
(2.6) yields the A Q for (2.1'), and using the H- inner product yields the A Q 
for (2.1), which is M^ 1 following the A Q for (2.1'). The Aq for (2.1’) is 
quite simple, and the reader might write it out. We will not, so that no one 
will think that we use it. We will point out that D(A Q ) requires both the 
geometric boundary conditions in V and the natural boundary condition w''(t,V) 
= 0; i.e. , zero moment on the right end. 


Ttemanic 6.J,. That the geometric boundary conditions 
w(t,0) = w'(t,0) = 0 


and 


( 6 . 12 ) 


w(t,J) = Wj^t) 


(6.13) 


are imposed in V but not in H — i.e., on the generalized displacement but not 
on the generalized velocity — is common in distributed models of flexible 
structures. The natural norm for expressing the kinetic energy of distributed 
components is the L^-norm, which cannot preserve constraints on sets of zero 
measure. Because the strain energy involves spacial derivatives, the stronger 
strain-energy norm can preserve the geometric boundary conditions (although, 
as for the boundary slope of an elastic plate, the V-norm may impose some of 
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these boundary conditions in an rather than a pointwise sense). The 
strain-energy norm is based on the material model of the distributed com- 
ponents of the system, and it should not be surprising that such a norm is 
required to connect the various structural components. Q 



6.2 The Optimal Control Problem 

We take Q = I in the performance index in (3.1). This means that the 
state weighting term <Qz,z>g is twice the total energy in the structure plus 
the square of the rigid-body rotation. Since there is one input, the control 
weighting R is a scalar. 

According to (3.33), the optimal control has the feedback form 
u(t) = - <f , x (t ) > v - <g , x( t ) >., 

(6.16) 

where x(t) has the form (6.1), and 

f = (a f ,0 f ,p f ) = R _1 n i B 0 e V, 

(6.17a) 
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(6.17b) 


g = ( WV = R 1 n i B 0 8H. 

Note that fi f = 0 f (|) is not used in the control law— recall (6.8) and (6.9). 
6.3 Approximation 


Our approximation of the distributed model of the structure is based on a 
finite element approximation of the beam which uses Hermite cubic splines as 
basis functions (I3X.W). These are the basis functions most commonly used in 
engineering finite element approximations of beams. The splines and their 
first derivatives are continuous at the nodes. Because the basis vectors e 

in the approximation scheme in Section 4 must be in the space V defined in ^ 
(6.10), we write them as 


5 1 = ( 1 , 0 , 0 ), 


(6.18a) 


B j - <°>0j.0j(l)), j = 2, 3 


» * ♦ t p 


(6.18b) 

where the dj's are the cubic splines. When we use % elements to approximate 
beam, there are 2n 0 linearly independent splines. Thus, with the rigid- 

body mode, the order of approximation is n = 2n + l 

e * 

For the numerical solution to the optimal control problan. we have only 

to plug into the formulas of Section 4. The matrices in (4.3) are calculated 
according to (4.3). < 4 . 3 , and (4.9). with B 0 glven by (s . 5) . 


K " U<e l' e J )1 ’ D “ * e 0 K n . H” - t<e 1 .e J > H ] 


(6.19) 


B 


'° - [10 0 ••• 01 - [<e 1> ^( l , 0 .0)) a ] . [<e (1 , 0<0)> j 

0 


( 6 . 20 ) 
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Note that the first row and column of K n are zero. The matrix!" in (4.32) is 
K n with 1 added to the first elanent. The matrices s" and B» are given by 
(4.5) and. since Q - I, the matrix Q" is the w" in (4.33). Kith these 
matrices, we solve the Rlccati sanation (4.40) and use (4.41) and (4.45) to 
compute the approximations to the functional gains, which are 


f n = Wn’^fn'Pfn* 


(6.21a) 


6 n “ (a gn'^gn' P gn > 


(6.21b) 


For convergence, we satisfy all the hypotheses of Theorem 5.9. In par- 
tlcular, since Q is the identity on E, it is ooerdve. Theorem 5.9 implies 
that the solutions to the finite dimensional Bicoatl sanations converge as in 
Theorem 5.1 and that the functional control gains converge as in theorem 5.10 

„,,„rW 6.2. It might appear that the hypotheses of Theorem 5.6 hold with n 0 , 
1, but not so. For J 1 2. e^ is orthogonal to e t in and V not in H. 

Recall (6.1>-<6.3>, (6.9>-(6.U> and (6.18). Q 
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6.4 Numerical Results 


Figures 5.2a and 6 . 2 b show the oomputed functional gain kernals g” and 
Ogn for the damping eoeffiolent c„ - 1<T 4 . the control weighting R . l and 

“e - 2. 3. 4. 5 and 8 beam elements. Table 6.2 lists the corresponding scalar 

components of the gains. For e 0 . lo ' 4 and R - 05. the convergence is slower. 

as discussed belt*. To sh« the complete story of convergence. Figures 6.3a 

and 6.3b and Table 6.3 shew the results for n - 2 3 4 5 a „ 

e 2 *4*5,8, and Figures 6 . 4 a 

and 6.4b and Table 6.4 shew the results for n e = 4,6,8,10. 

We have plotted d fn because the second derivative appears in the 
strain-energy inner product in < 6 . 8 , and < 6 . 9 , and d fn converges in 
Note that, since the Hermlte cubic splines have discontinuous second deriva- 
tives at toe uodes, toe approximations to are discontinuous at the nodes. 

Although H -convergence guarantees onlv I ' • 

8 S onl y ^-convergence for 0 fn , it can be 

shown that g” converges uniformly on [O.jf] for this problem. 

UtleS Mlt »*. *> -*•*«• the fact that it does not appear in the 
feedback law and the fact that the convergence of Pf# is not an independent 

Piece of information about the convergence of toe control gains; since d fn <0, 

Pf n ' toe convergence of dy n implies the convergence of 

On the other hand, although 8 = * <0\ « 

sn ^gn W for each n, the H-norm convergence 

Of g„ does not enforce this condition in the limit, as the V-norm convergence 

of f„ enforces S f - d f (|, . Hence, as far as we can tell from our results in 

Sections 3.5, Sfn is an independent indicator of toe oonvergenoe of toe con- 
trol gains, as well as being used in the control la, in <6.16,. However, toe 
behavior of dgn in Figures 6 . 2 b. 6.3b and 6 . 4 b suggests that g„ converges in 
V. Stronger results on the continuity of 0g and the convergence cf d (jf, 
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should follow frcm a theory stating that, because the ope^loop semigroup 
generator 1 is analytic, the solution to the infinite dimensional Ricoati 

equation naps all of E into DU*). £<*> «»" ** * ”” ^ 

Flg ure 6. 2a also suggests such a theorem, but we have not proved it. 

With the state weighting Q fixed, the two factors that determine the rate 
of convergence are c 0 and R. Mthou*. we have used splines to approximate the 
beam, the relation between the convergence rate and c„ and R probably can be 
interpreted best in terms of the number of .tural modes of the structure that 
the optimal infinite dimensional controller really controls. Strictly speak 
!„g. the controller controls all modes, but the functional gains lie essen- 

tially in the span of some finite n-ber of modes. This would be the n^ber 

- Mins if we used the natural modes as 

of modes required for convergence of the gain 

the basis vectors in the approximation. Tbe rest of the modes are practically 

(but not exactly) orthogonal to the functional gains, so that the optimal 

a Tn eeneral, the lighter the damping, 

feedback law essentially ignores than. In gener 

a „ that will be controlled for given Q and R; the cheaper the con- 
th6 more modes t-het wii-L 

, ... Blll controlled for given Q and c Q . The question 

trol, the more modes that will oe conuun 

of the convergence of the finite eluent ap^oximation to the functional gains 

becomes then a question of hew «a W modes the optimal control law really wants 

im»«xxni-<» it takes to approximate those modes, 
and how many elements it taxes w 

Numerical experience with optimal control of flexible structures has 
shown this modal interpretation of the convergence of the approximating con- 
trol laws to to very useful, and that it is difficult to Improve upon the 
natural modes as basis vectors for the approximation scheme (see IG5» «»- 
ever, whether the natural modes are always or almost always the best basis 
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vectors is an open question. We use the cubic splines here to demonstrate 
that a standard finite element approximation works quite well. Also, to use 
the natural modes as basis vectors here, we first would have to compute them 
using a finite element approximation -- as in most real problems — and we do 
not know in advance which or how many modes are needed. On the other hand, if 
the most important natural modes are determined from experiment, then modal 
approximation should be best. 
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Figure 6 ,2a. Functional Control Gain Component i {n ' ' 
Damping coefficient c Q = 10 -4 ; control weighting R = 1 
number of elements n 0 *= 2, 3, 4 , 5, 8 


i^Too FOo 3b. oo ibToo sb.oo so.oo 70.00 io.oo *0.00 io< 

Figure 6. 2Ji. Functional Control Garin Component 

Damping coefficient c Q = 10 -4 ; control weighting R - 1 

number elements n g * 2 # 3, 4, 5, 8 






cn 

to 



number of elements n fl = 2, 3. 4, 5, 8 


155.00 



number of elements n g = 2. 3. 4. 5. 8 





o> 

<n 



nuaber of lento . 4 , 6 . J( , 


Flg„r. ""‘ * ” ” .. .. .. . 

-■-i. Functional Control Gun Component t 

°“ Pln8 C ° , " 1Sl “ t «. * 10-; control w.^tlns “ 
nun bar of .lento , 4 . 



n 

e 

°fn 

a_„ 

gn 

^gn 

2 

1.000 

.2141 

-22.459 

3 

1.000 

.23 96 

-25.221 

4 

1.000 

.2496 

-26.331 

5 

1.000 

.2534 

-26.786 

8 

1.000 

.2561 

-27 .041 

Table 6.2 

. Scalar 

Components 

of Functional Control Gains 

Damping coefficient 

O 

o 

II 

M 

O 

1 

control weighting R = l 

number of 

elements 

n e = 2 ’ 3 ' 

00 

■n 

«n 

n e 

a fn 

a gn 

^gn 

2 

4.4721 

1.0136 

-108.27 

3 

4.4721 

1.1770 

-126.40 

4 

4.4721 

1.2440 

-133.87 

5 

4.4721 

1.2781 

-137.57 

8 

4.4721 

1.3106 

-141.15 

Table 6.3. 

Scalar Components of Functional Control Gains 

Damping coefficient 

• «* 

1 

o 

rH 

II 

o 

o 

control weighting R = .05 

number of elements n 

e 2,3,4 

. 5, 8 

n 

e 

°fn 

“gn 

^gn 


4.4721 1.2440 -133.87 

4.4721 1.2973 -139.69 

4.4721 1.3106 -141.15 

4.4721 1.3141 -141.54 


Iafel§ 6.4. Scalar Components of Functional Control Gains 
tempins coefficient c 0 , 10 -t. uelghtlng R , 05 

number of elements n g = 4, 6, 8, 10 
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Figures 6.5a and 6.5b and Table 6.5 represent attempts to compute an 
optimal control law for the structure when R = .05 but c Q = 0. Since Q is the 
identity operator in E and hence coercive. Theorem 5.4 says that no optimal 
control law exists and that the norm of the solution to the finite dimensional 
Riccati equation grows without bound as the number of elements increases. 

This is reflected in the nonconvergence of a g n » 0gn an< * ^gn' fn 

verges and the convergence of 0 fn is unclear. 

in applications where the structural damping is not known, except that it 
is very light, it is tempting and not uncommon engineering practice to assume 
aero damping in the design of a control law for the first few modes, while 
trusting whatever damping is in the higher modes to take care of them. How- 
ever. if high performance requirements (large Q) or coupling between modes in 
the closed-loop system necessitate a control law based on a more accurate 
approximation of the structure. Theorem, 5 .4 and the current example warn that 
the higher-order control laws are likely meaningless and rather strange if no 

damping is modeled. 

W. should note that we have seen similar problems [G91 where n„ remains 
bounded and the gains oonverge for zero damping but finite-rank Q. In such 
cases. Theorem 5 .3 says that an optimal control law exists for the distributed 
model of the structure and that the finite dimensional control laws converge 
to an optimal infinite dimensional control law. Also. Balakrishnan IB11 has 
shown that an infinite dimensional optimal control law exists for no damping 

when Q = BB*. 
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n e 

°fn 

V 

Pgn 

2 

4.4721 

1.0516 

-112 .23 

3 

4.4721 

1.3061 

-140.18 

4 

4.4721 

1.4758 

-159.11 

5 

4.4721 

1.5996 

-172.64 

8 

4.4721 

1.8407 

-199.39 


Table 6.5. Scalar Components of Functional Control Gains 
Zero damping; control weighting R = .05 
number of elements n g = 2, 3, 4, 5, 8 
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7 . The Optimal Infinite Dimensional Estimato r, Compensator and Cl££ed-l<2fi£ 
System 

As in Sections 3 and 5, we will state seme initial definitions and 
results for an arbitrary linear control system on a Hilbert space in Subsec- 
tion 7.1, and then discuss implications for flexible- structure control in Sub- 
section 7.2. 

7.1 The Generic Problem 

Let A, T(t) and B be as in Subsection 3.1, with E an arbitrary real Hil- 
bert space. The differential equation corresponding to (3.2) is, of course, 

z(t) = Az(t) + Bu(t), t > 0. 

(7 .1) 

We assume that we have a p-dimensional measurement vector y(t) given by 

y(t) = C Q u(t) + Cz ( t ) , 2) 

where C Q e L(R m , R P ) and C e L(E,R P ) for some positive integer p. 

Definition 1.1. For any F e L(R P ,E), the system 

• ^ 

z(t) — Az (t) + Bu(t) + F [y(t ) — C n u(t)-Cz(t)l , t > 0, 

u (7.3] 

will be called an observer , estimator (we use the terms interchangeably) , for 

A A 

the system (7.1)-(7.2). Let S(t) be the semigroup generated by A-FC. The 

A 

observer in (7.3) is strongly (uniformly exponentially) stable if S(t) is 
strongly (uniformly exponentially) stable. □ 

To justify this definition, we write 
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e(t) 


z(t) - 2(t) 


(7.4) 


and, with (7.1)-(7.3), obtain 

e(t) = 3(t)e(0), t 2 0. 


(7.5) 


Of course, an observer, or estimator, is necessary because the full state 
z(t) will not be available for direct feedback, and the feedback control must 
be based on an estimate of z(t). When, as in this paper, the desired control 
law has the form 


u(t ) = -Fz ( t ) 


(7.6) 


for some F e L(E,R m ), the observer in (7.3) can be used to construct z(t) from 
the measurement in (7.2) and then the control law in (7.6) can be applied to 
z(t). The control applied to the system is then 


u ( t ) = -Fz ( t ) , 


and the resulting closed-loop system is 



t 2 0, 


where S^.^tt) is the semigroup generated on E x E by the operator 


A -BF ] 

$C [A-BF-$C]J' D(A -'»> = D(A)XD(A). 


(7.7) 


(7.8) 


(7.9) 


With the estimator error e(t) defined by (7.4), it is easy to show that 

(7.8) is equivalent to (7.5) and 
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z(t) 


(A-BF)z(t) + BFe(t), t > 0, 


(7.10) 


where (A-BF) generates a semigroup S(t) on E. Also, it is easy to prove the 
following. 

Theorem 7.2. Suppose that there exist positive constants M^, Mjj, a 2 and Oj 
such that 

t 

I IS(t) | | < e 


-o 2 t 


||S(t)|| < V^e , t 2 0 . 


(7.11) 


Then 


, for each real a, < minUj, o 2 ) , there exists a constant M3 such that 


II Wt>M ^3 e 


-a 3 t 


, t 2 «». 


(7.12) 


Also, 


a(k,J = o(A-BF) U ff(A-FC), 


(7.13) 


where o(A ) is the spectrum of A^,^. P 


The observer in (7.3) and the control 1?« in (7.7) constitute a compensa 
tor for the oontrol system in (7.1) and (7.2'. The transfer function of this 
compensator is 

fi(s) = -F (si- [A-BF + F(C F-C)]) _1 f , 

u (7*14; 

which is an m X P matrix function of the complex variable s. When E has 
infinite dimension, the compensator transfer function is irrational, except in 
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degenerate, usually unimportant cases. 


The foregoing definitions of this section and Theorem 7.1 are straight- 
forward generalizations to infinite dimensions of observer-controller results 
in finite dimensions. Balas [B3] and Schumacher [S2] have used similar exten- 
sions. 


Now suppose that F is chosen as 

f ■ 

(7.15) 

A 

where II c L(E,E) is the minimal nonnegative selfadjoint solution to the Ric- 
cati equation 

AH + HA* - n C* R _1 C n + $ - 0 . 

(7.16) 

A /V 

with Q e L(E, E) nonnegative and selfadjoint and R e L(R p ,R p ) symmetric and 

positive definite. Theorem 3.3 (with A, B, Q, R, Hand S(t) replaced by 

* * A A A A * a 

A ,C , Q, R, H, and S (t)) gives sufficient conditions for ft to exist and for 

A £ A 

the semigroup S (t) — and equivalently its adjoint, the S(t) generated by 

d * 

A ll c R c — to be uniformly exponentially stable. 


Definition When the control gain operator is 

F = R _ 1 B*n , 


(7.17) 


with II the solution to the Riccati equation (3.3) , and the observer gain 
operator is given by (7.15) and (7.16), we will call the compensator consist- 
ing of the observer in (7.3) and the control law in (7.7) the infinite 

dlffien3l,Qn3l compensator, and (7.8) the optimal closed-loop system . □ 
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Optimal Infinite Dimensional Compensator 


Figure 7 .1 . Optimal Closed-loop Sj stem 


Remark 7.4. The infinite dimensional observer defined by (7.3), (7.15) and 
(7.16) is the optimal estimator for the stochastic version of (7.1) and (7.2) 
when (7.1) is disturbed by a stationary gaossi an white noise process with zero 
mean and covariance operator Q and the measurement in (7.2) is contaminated by 
similar noise with covariance R. For infinite dimensional stochastic estima- 
tion and control, see [Bl, C4] . When the state weighting operator Q in (3.1) 
is trace class, the optimal infinite dimensional compensator minimizes the 
time-average of the expected steady-state value of the integrand in (3.1). 
Existing theory for stochastic control of infinite dimensional systems 
requires trace-class Q, but we have a well defined compensator for any bounded 
nonnegative selfadjoint Q and $, as long es the solutions to the Riccati 
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equations exist. As the next two sections show (without assuming trace-class 

A 

Q) , the infinite dimensional compensator is the limit of a sequence of finite 
dimensional compensators, each of which can be interpreted as an optimal LQG 
compensator for a finite dimensional model of the structure. Therefore, we do 
not require trace-class Q in our definition of the optimal compensator, even 
though this compensator solves a precise optimization problem only when Q is 
trace class. 

This paper is concerned primarily with how the finite dimensional compen- 
sators converge to the infinite dimensional compensator, and the a na lysis of 
this convergence requires only the theory of infinite dimensional Riccati 
equations for deterministic optimal control problems and the corresponding 
approximation theory. While the stochastic interpretation of the infinite 
dimensional compensator and, in Section 8.2, of the finite dimensional estima- 
tors should be motivational, nothing in the rest of the paper depends on a 
stochastic formulation. We assume that the operators Q, R, Q and R are deter- 
mined by some design criteria. In many engineering applications, determinis- 
tic criteria such as the stability margin and robustness of the closed-loop 
system, rather than a stochastic performance index and an assumed noise model, 
govern the choice of Q, R, Q and *.Q 

7*2 Application to Structures 

For the rest of the paper, E = V X H as in Section 2, and A and B are the 
operators defined there. 

The measurement operator C in (7.2) now must have the form 
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c = [C 1 C 2 ] 


(7.18) 
th 


where C ± e L(V.R P ) and C 2 e L(H,R P ). Hence, if we denote by (C(x.x)) i the i 
component of the ^-vector C(x,x> , for (x,x) e E, then there must exist c u e V 
and 8 H such that 


(CCx.x)^ = <C ljL ,x> v + <C 2 i' 3> H' 1 " 1, **'* p * 


(7.19) 


Also, the estimator gain operator F is givnn by 


a P a A 

Fy = I (f,,g 1 )y i 

i=l 1 


( 7 . 20 ) 


for y = [y 2 y 2 ... y p 
are elements of V and H, respectively 


A A 

] T 6 R p , where the functional egtAmatQr gaAliS fj_ and gj^ 


For the optimal estimator gains, we cm partition II as 


«£ ft 

n» ni 

A* A 

Hi n. 


( 7 . 21 ) 


and use (7.15) and (7.19) to get 


A 

f. 


I (R + • 


(7.22a) 


A 

«i 


I (R -1 ) ( ^i c i j + D2 c 2j ) ' 1 = 1,2 


(7.22b) 


Now let us partition Q as in (4.34): 
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( 7 . 23 ) 


$ = 


rA A 

Q 0 Q 1 

A* 4 


In the optimal control problem, we almost always have a nonzero Q Q because 

this operator penalizes the generalized displacement. For the results in this 

A 

paper, Q Q can be nonzero in the observer problem, and, as in the control prob- 
lem, some of the strongest convergence results for finite dimensional approxi- 
mations can be proved only for coercive $. However, if the observer is to be 
thought of as an optimal filter, then Q should be the covariance operator of 
the noise that disturbs (2.1). l n this case, Q Q = o and = 0. 


78 



8. Approximation of the Infinite Dimensional Estimator 


8.1 The Approximating Finite Dimensional Estimators 

Here, the scheme for the approximation of the flexible structure is that 
in Section 4. We will construct on the sutspace an estimator that approxi- 
mates the optimal infinite dimensional estimator of Section 7, and this esti- 
mator will produce an n th -order estimate z r = of infinite dimen- 

sional state vector z = (x,x). In Section 9, the the n -order compensator 
that results from applying the n th approxin ation to the optimal control law 
(in Section 4.2) to z r will approximate the optimal infinite dimensional com- 
pensator of Section 7. 

H^afcheais 8.1. There exist a sequence e L(E n ,R P > such that 

H C n P En- C " -*• ° aS (8.1) 

and a sequence Q n e L(E n > , Q R = Q n 2. 0* inch that 

°n P En S strongly as n -» - -O (8 . 2 > 

Hypothesis 8.2. For each n, the system (A^.CM is stabilizable. In particu- 
lar, any unstable modes of the system (C^, are observable. D 


The n th observer, or n u “ estimator , 1.3 


th 


S n = A n*n + B n u + F n (y-C 0 u_C n^n ) 


(8.3) 


where the estimator gain F r is 
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(8.4) 


A /V *A_1 

F n ' n n C n R 1 


and is the nonnegative selfadjoint solution to the Riccati operator equa- 
tion 


a A + - V?" lc A * K ■ o- 

Hypothesis 8.2 implies that such a solution exists and is unique. 


(8.5) 


This representation of the n^* 1 estimator as a system on E , with the 
estimator gain determined by the solution to a Riccati operator equation, is 
necessary for showing how the sequence of finite dimensional estimators 
approximate the infinite dimensional estimator. However, on-line computations 
will be based on the equivalent differential equation 

H = A n $ + B n u + F^y-C^u-C 1 ^) 

( 8 . 6 ) 

where ij(t) 8 R^ n , A n and B n are the matrix representations of the operators A r 
and B n , as in Section 4, and C n is the matrix representation of C . 

A_ 

The 2nxp gain matrix F is 

F n = fl” VT n (C n ) T R -1 , 

(8.7) 

- A 

where W is the 2n X 2n grammian matrix in (4.33) and U* 1 satisfies 

A n n° + ft 1 w _n (A n ) T w n + n * 1 w _n (c n ) T R _ 1 c n n n + o n = o, 

( 8 . 8 ) 

with Q n the matrix representation of The relationship between z = 

jA A 

t* n .x n ) and tj is, of course. 
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(8.9) 


x n (t) = I « i (t)e i 
n i=l 


and 


Since the matrix representations of A n and A R are A n and W n (A ) W , respec 
tively, and the matrix representation of C* is W n (C n ) T , (8.7) is the matrix 

representation of (8.4), and the 2nX2n Riccati matrix equation (8.8) is the 

A A 

matrix representation of (8.5), with IP 1 the matrix representation of D^. 
(Recall that W -n is the inverse of W n .) 


As in the control problem, we do not solve the matrix representation of 
the n th Riccati operator equation directly because the matrix representation 
of a self ad joint operator in general is not symmetric. In the duality between 
the optimal control and estimator problems, (8.5) and (8.8) correspond to 
(4.26) and (4.38). respectively. In (4.39), we defined the symmetric matrix 
5“ = W n rf 1 and then obtained the Riccati equation (4.40) to solve for 5” . 

We proceed in a similar fashion here, but with an interesting difference. 


Since II and <5 are nonnegative selfadjoint operators on E 
n n A 

Q n are their matrix representations, the matrices W n ft 3 and 
tive and symmetric. Hence, the matrices 


and If 1 and 
n 

are nonnega- 



( 8 . 11 ) 


and 
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A 

Q 


n 


= qV" 


are nonnegative and symmetric. 
(8.8) yields 


( 8 . 12 ) 

Substituting (8.11) and (8.12) into (8.7) and 


A 

F 


n 


= fr(C n ) T 



and 


(8.13) 


A n fl“ + n"(A") T - &'<C n > T R- 1 C n l5" * Q n - 0. 

(8.14) 

the Riccati matrix equation to be solved numerically in the n th approximation 
to the infinite dimensional estimator. In view of the relationship between 
(8.5) and (8.8) and the relationship between (8.8) and (8.14), we see that 
Hypothesis 8.1 guarantees the existence of a unique nonnegative symmetric 
solution to (8.14). 


To see the relationship between the matrices in (8.14) and the operators 
in (8.5) more clearly and the difference between the current approximation 


scheme and that used in Section 4.2 for the control problem — suppose that we 

^ A 

take Q n — Pg^Q|g^. Let Q be defined as in (4.36) and (4.37) with Qq, , and 
Qj replaced by Q Q , , and • Then 

*n . W -n?n w - ni 

(8.15) 


For example, if Q in the control problem and Q in the estimator problem are 
both equal to the identity, then the Q n in (4.35) - (4.42) is W n and 

/V 

A n — j-j 

Q = W . This may seem suspicious, but Subsection 8.2 should demonstrate 
that we are solving the appropriate estimator problem here. 


The only thing missing now for numerical implementation of the n th 
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estimator, or observer, is to give C n , the matrix representation of C n . expli- 
citly. We write 


r,n 


= [C? C?1 


(8.16) 


where the p x n matrices C» and c" are, respectively. the matrix representa- 


tions of the operators C ± and Cj in (7.18). We can cover virtually all appli- 
cations hy assuming C n = C| En> in which case the 1 th colvmn cf cj is the p- 
vector equal to C^. and the i th column of is the p-vector equal to C^. 

*s mm ims iis usaBlsis asi si sasall:m Isl naasEisal inplenentaU ai) si 

ite n th state estimator ; For online computation, the n th estimator, or 
observer, is (8.6). the gain matrix?” is given by (8.13) and the solution to 

<9f 

. . / q . . \ Th ft matrices and C are defined as 

the Riccati matrix equation (8.14). Ihe matrices u 

above. 

8.2 Stochastic Interpretation of the Approximating Estimators 

As we have said, our approximation theory for the optimal estimator is 
based on approximation of the infinite dimensional Hiecati equation, whose 
structure is the same for both control and estimator problems, and the sto- 
chastic properties of the optimal estimator problem never enter our approxima- 
tion theory. Furthermore, using only the deterministic setting above, we will 
proceed, subsequently, to analyze the finite dimensional estimators and the 
compensators based upon them. Ncntheless. we should consider momentarily the 
sequence of finite dimensional stochastic estimation problems whose solutions 
are given by the equations of the preceding subsection. 

First, recall how the covariance operator of a Hilbert space-valued ran- 
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dom variable is defined. The covariance operator of an E-valued random vari- 
able (i> is the operator Q for which 


expected value Kz,u>> E <S,o» E } = <0z.%> E , z, z e E. 

(8.17) 

(See [Bl, C4] .) 

With F n given by (8.4) and (8.5), (8.3) is the Kalman-Bucy filter for the 
system 


(8.18) 

y * c 0 u + c „ z n + "O' 

(8.19) 

where w R (t) is an E n ~valued white noise process with covariance operator Q n 
and <*>Q(t) is an R^-valued white noise process with covariance operator 
(matrix) R. Next, careful 1 inspection will show that the filter defined by 

(8.6), (8.13) and (8.14) is the matrix representation of the filter defined by 
(8.3) , (8.4) and (8.5) . 


With z n and n related as in (4.1) and (4.4), (8.18) and (8.19) are 
equivalent to the system 


n = A n n + B n u + (/, 

y = c 0 u + c n n + u) 0 , 

where (/( t) is the R^ n -valued noise process related to o> (t) 


( 8 . 20 ) 


( 8 . 21 ) 


n (t > by 
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( 8 . 22 ) 


<i) (t) = 21 ((A(t)e^> • 

n i=l 

Certainly, a Kalman-Bucy filter for (8.20) and (8.21) has the form (8.6) with 
the filter gain given by (8.13) and (8.14). This particular filter is the 
matrix representation of the filter defined by (8.3), (8.4) and (8.5) if and 
only if the matrix Q n defined by (8.12) is the covariance of the process (/(t). 
Since Q n is the matrix representation of . straightforward calculation using 
(8.12) and (8.17) shows that the* n in (8.12) is indeed the correct covariance 

matrix. 

Of course, if « n (t) and (/(t) represent a physical disturbance to the 
structure, then % <t> must have the form (0,m< 2) (t)) and the first n elements 
of (/(t) must be zero, but this is not necessary for our analysis. 

The finite dimensional observers can be interpreted now as a sequence of 
filters designed for the sequence of finite dimensional approximations to the 
flexible structure, with the n th approximate system disturbed by the noise 
process o> n (t), whose covariance operator is Q n « By hypothesis 8.1, these 
covariance operators converge to the operator of Section 7. If we have a 
reliable model of a stationary, zero-mean gaussian disturbance for the struc- 
ture, then we can take the covariance operator for this disturbance to be Q 
and think of the infinite dimensional observer as the optimal estimator. But, 
again, this interpretation is not necessary for the rest of our analysis. 

8.3 The Approximating Functional Estimator Gains 

The n th estimator gain operator in (8.4) has the same form as the infin- 
ite dimensional estimator gain in (7.15) and (7.20). We have 
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A P A - 

V = £ j 


i=l 


in’ g in ,y i 


(8.23) 


for y [y 1 y 2 ••• y p J T e R p , where the functional estimator gains £ and 
g in are elements of V R = H n> The matrix F n in (8.7) and (8.13) is the matrix 

A 

representation of F n# which means that, if we write 


F 


f l f 2 ... 


P P 

g l 

LP 1 P 


Si &2 

X O X • « • 


Sr 


f g 

where the columns P i , p 1 e R n , then 


(8.24) 


in 


n f ± 

Z \ 
j=l J J 


i 1. • • • f Pi 


(8.25a) 


a n g ± 

g in ^ i = 1# • • . »p. 

j=l J J ' 


(8.25b) 

For convergence analysis, it is useful to note that f ±n and g in are given 
also by equations corresponding to (7.22). With the measurement operator C 
written as in (7.19) and C n = c| E , we have 


,*-l 


in 


J ? 1 (s hon P Vn c ij + 


(8.26a) 


where 


j Z l <R 1> ij ( n in P Vn c ij + ^n^n^' 


in 


(8.26b) 


A 

“n = 


A 

^n 

a* 

^n 


A 

n 


In 

in 


(8.27) 
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8.4 Convergence 


Now we will indicate the sense in which the finite dimensional 
estimators/observers approximate the infinite dimensional estimator in Section 
7. As we have said, implementation of the n th estimator is based on (8.6), 

(8.13) and (8.14), but convergence analysis is based on the equivalent system 
(8.3), (8.4) and (8.5). The question then is how the observer in (8.3), with 
gain given by (8.4) and (8.5), converges to the observer in (7.3) with gain 
given by (7.15) and (7.16). 

Recall Hypothesis 8.1, and recall fron Section 4 that the approximations 
to both the open- loop semigroup and its adjoint converge strongly. Also, 
recall that Hypothesis 8.2 guarantees a unique nonnegative selfadjoint solu- 
tion to the Riccati equation (8.5) for each n. Replacing A r and B n with A R 

and C* in Theorems 5.1 and 5.3, we obtain 
n 

Theorem 8.3. i) If II EL 1 1 is bounded uniformly in n, then the Riccati 

“ A A 

algebraic equation (7.16) has a nonnegative selfadjoint solution II and HoPeh 
converges weakly to ft. ii) If there exist positive constants M and p, 
independent of n, such that 

llexpItAn-f^V^nH’ 11 i Me' 111 • *20, (8>28) 

A A 

then || || is bounded uniformly in n, U n P En converges strongly to II and 

exp( lA n ~ D n C*R - ^C n ]t )P En converges strongly to§(t), the semigroup generated 
by A-ftc*R~*C, the convergence uniform in t 2. 0 . iii) If is bounded away 
from zero uniformly in n, then llfi^ll being bounded uniformly in n guarantees 
the existence of positive constants M and P for which (8.28) holds for all n. □ 
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The proof of the following theorem is practically identical to that of 


Theorem 5 .4 . 


Iil£2Eea 8.4. If Q is E-coercive and d Q = 0 , then there is no nonnegative sel- 
fadjoint solution of the Riccati operator equation (7.1 6), and 

A 

IlDnll 00 as n — > °°. □ 

(8.29) 


Our purpose for bothering to state this obvious dual result is to point out 
the following question. Can Theorem 8.4 be modified to include the case where 
Q has the form (7.23) with Q Q = 0 , \ = 0 and Q 2 coercive on H? 


Next, we have the dual result to Theorem 5.6: 


XteP . rew 8.5, Suppose that has an invariant subspace Vq which is also invari- 
ant under the damping map that E Q = V Q X V Q is an observable subspace, 

and that the restrictions of A Q and d 0 (*,O to are both H-coercive. Also, 
suppose that V Q has finite dimension n Q and that, for each n 2 in the 
approximation scheme, the first Hq e i 's span V Q and the rest are orthogonal to 
V Q in both V and H. 

i) Then (7.16) has nonnegative solution H, and || TT^ || is bounded uniformly 
in n, so that ft^P^ converges to ft weakly. 

ii) If E Q and E^* (the E-orthogonal complement of E^) are invariant under Q, 
and if the E Q -part of the astern (A,Q) is controllable, then the hypothesis of 
Theorem 8.3 ii) holds. 


JiEflaf. The proof is practically identical to that of Theorem 5.6 with B 
replaced by C . For ii) , note that, when we partition A and ^ as in (3.16), 
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the finite dimensional system is controllable if and only if the 

system ^i’Ajj) ls observable. □ 

Remark 8.£. Remarks 5.7 and 5.8 pertain to Theorem 8.5 as well as to Theorem 
5.6; i.e. , in most applications the theorem requires either that both A Q and 
d Q be coercive (so that V Q = {0}) or that the natural mode shapes be the basis 
vectors and the damping not couple the natural modes. It seems unlikely that 
a finite number of observable rigid-body modes could change the nature of the 
convergence, but they greatly complicate the proofs. For applications where 
both rigid-body displacement and rigid-body velocity are measured, a result 
analogous to Theorem 5.9 can be obtained, but we will not bother here because 
it adds no significant insight and we cannot use it in the example in Sections 

6 and 10. Also, see Remark 10.1. □ 

Theorem 8.7. If converges strongly toll, then 

I l f in " h “*° (8.30a) 

I lg in - gjj l H -* °* as m * 1 

where f^ and g^ are the functional estimator gains 
are the approximating functional gains in (8.25). 

Proof . The result follows from (7.22) and (8.26). 


(8.30b) 

in (7.20) and £ ±n and g in 
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9. The Finite Dimensional Compensators and Realizable Closed-loop Systems 
9.1 Closing the Loop 

The n th compensator consists of the n th approximation to the optimal con- 
trol law in Section 4. applied to the output of the n th estimator/observer in 
Section 8; i.e*, the feedback control 

u = — f £ 

n r n z n 

(9*1) 

where 

F n “ «“ lB X 

(9.2) 

(recall (4.25)) and z n (t) is the solution to (8.3). Equivalently, this com- 
pensator can be written as 


u„ - - 


where 


(9.3) 


F n = R 1 B nT ff 1 

(9.4) 

(recall (4.43)) and the 2n-vector ^(t) is the solution to (8.6)). On-line com- 
putations will be based upon the latter representation, and the block diagram 

in Figure 9.1 shows the realizable closed-loop system that results from the 
bh 

n compensator. We will refer to this system as the i£ h closed-loop system . 
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Control System 



n th com pe nsa tor 

Figure 9.1 Closed-loop System 


This closed-loop system is equivalent to 


A 

z 


A 

Z 


where the operator 


k= Da *- n) ■ D<A,XE - 


generates the closed-loop semigroup S w , n (t) on EX E n . The closed-loop 
response produced by the n th compensator — i.e. , the response of the i 


(9.5) 

( 9 . 6 ) 

,th 


closed— loop system 


can be written then as 




Note that k <oi ^ haa oompaot resolvent if and only if A does. 

9.2 Convergence of the Closed-loop Systems 

Now we will consider the sense in which the n th closed-loop system 
approximates the optimal closed-loop system in Section 7 (Definition 7 . 3 ). 
Recall from Sections 4.1 and 8.1 how the approximating open-loop semigroups 

T n (*) and their adjoints converge strongly and how the input operators B n , the 
measurement operators C r and their respective adjoints converge in norm. Sec- 
tions 5 and 8 have given sufficient conditions for the approximating control 
and estimator gains to oonverge to the gains for the optimal infinite dimen- 
sional compensator. In this section, we will assume 

Hypothesis 9 . 1 . As n — » », 

1 ,F n p En - FH-* °> 

(9.8) 

NF n -F||-»0. □ 

(9.9) 

Remark 9.2. Of oourse, we are interested primarily in the case where the 

A 

gains F and F are the optimal LQG gains in (7.15) and (7.17) and F q and £ ape 
the corresponding approximations in Sections 4 and 8 (i.e. ,( 9 . 2 ) and ( 8 . 4 )). 
However, for the analysis of this section, we need only Hypothesis 9.1 for 
some F e L(E,R m ), p e L(R P ,E) and approximating sequences F n and F n . Any such 
gain operators will yield olosed-loop semigroup generators A^,^ i n ( 7 . 9 ) a nd 
A .. n in (9.6). □ 
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We denote the projection of E X E onto E X E n by P EEn » 

Theorem 9.3. For t 20, s «» n (t)P EEn converges strongly to S^.^t). and the 
convergence is uniform in t for t in bounded intervals. 

Proof . This follows from the strong convergence of the open-loop semigroups 
and the uniform norm convergence of the control and estimator gains. □ 

We should expect at least Theorem 9.3, but we need more. We should 

require, for example, that if S(t) is uniformly exponentially stable, then 

S , (t) must be also for n sufficiently large. Although nunerical results for 
00 n 

nunerous examples with various kinds of damping and approximations suggest 
that this is usually true, we have been unable to prove it in general. We do 
have the result for the following important case. 

Suppose that the basis vectors 6j of the approximation scheme are the 
natural modes of undamped free vibration and that the structural damping does 
not couple the modes. Then, for each each n, E n and E^ reduce the open- loop 
semigroup T(t) and its generator A. For this case, we can extend A^, n to D(A)X 

D(A) as 

_ r A -BF P Pn l 

X = U /r n En 

®’n |F n C n Comp J (9.10) 

where 

^nComp = [ V B n F n" F n C n lP En + A, D(A>n £ (? n) 

Note that Ejf"is the span of the modes not represented in the n th compensator. 
The operator A^, generates a semigroup S <D , n (t) on E X E. X E fl and 
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(OJXEjJ’ reduce S^.^t), and the restriction of s’, (t) to EXE is S , (t). 

n n ® n 

A/ 

Hence s m * n (t) is uniformly exponentially stable if and only if both # n ( t ) 

and the part of the open-loop system on EjJ* are uniformly exponentially stable. 

T freQrem £*4. i) Suppose that the basis vectors of the approximation scheme 
are the natural modes of undamped free vibration and that the structural damp- 
ing does not couple the modes. Then S . (t) converges in norm to S , (t), 
uniformly in bounded t-intervals. 


ii) If, additionally, S^.^t) is uniformly exponentially stable, then S,,, (t) 
is uniformly exponentially stable for n sufficiently large. 

Proof . From (9.6), (9.10) and (9.11), we have 


where 



0 

A /\ 

[F-F n JC 


B1F n p E,f FI | 
A n J 


(9.12) 


4 n • (B „ F n P En- BF) + (F n C n P En- FC > ‘ 

(9.13) 

Therefore, I l A a>'a» -A ®' n ll — > 0 as n — >®, and the theorem follows. □ 

This paper emphasizes using the convergence of the approximating control 
and estimator gain operators F n and F n , and the convergence of the functional 
gains that can be used to represent these operators, to determine the finite 
dimensional compensator that will produce essentially optimal closed-loop per- 
formance. However, close examination of the right sides of (9.12) and (9.13) 
reveals another important convergence question. While the gain convergence in 
(9.8) and (9.9) drives the off— diagonal blocks in (9.12) to zero, the norm 
convergence of the approximating input and output operators also is essential 

y-. 
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in killing A 


Expanding the two terms in this block yields 


n* 

B F Pn - BF = B (F P p -P ) + (B -B)F, 

n n En n n En n 


(9.14) 


A A 

F C Pp, - FC 
n n En 


A A 

(F -F)C Pp 
n n En 


F(c n p En- C) 


(9.15) 


The second term on the right side of each of' these equations represents, 
respectively, control and observation spillover, which has been studied exten- 
sively by Balas. Together, the control spillover and observation spillover 
couple the modes modelled in the compensator with the modes not modelled in 
the compensator. The spillover must go to zero — as it does when B r and C R 

A/ 

converge — for m n to go to zero* 


We should ask then whether there exists a correlation between the con- 
vergence of F and F and the elimination of spillover. The answer is yes if 
n n 

no modes lie in the null space of the state weighting operator Q in the per- 
formance index and if the assumed process noise, whose covariance operator is 
Q, excites all modes, but this correlation is difficult to quantify. As we 
discussed in Section 6.4, the two main factors that determine the convergence 
rates of the gains are the Q-to-R ratio and the damping, neither of which 

affects the convergence of B and C„. On the other hand, when either factor 

n n 

(small Q/R or large damping) causes the gains to converge fast, it generally 

A 

also causes the magnitude of F and F to be relatively small, thereby reducing 
the magnitude of the spillover terms in (9.14) and (9.15). Also, as n 
increases, the increasing frequencies of the truncated modes usually reduce 
the coupling effect of spillover. This is well known, although it cannot be 
seen from the equations here. In examples that we have worked, we have found 

that when n is large enough to produce convergence of the control and 
estimator gains, the effect of any remaining spillover is negligible. But 

this may not always be true, and spillover should be remembered. 
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9.3 Convergence of the Compensator Transfer Functions 


The transfer function of the n™ 1 compensator (shown in the bottom block 
of Figure 9.1) is 


*„<=> 


-V s1 - [ V B n F A< c 


O'n-V" 


- 1 * 


n’ 


(9.16) 


which is an m X p matrix function of the complex variable s for each n, as is 
the similar transfer function S(s) in (7.14) for the infinite dimensional com- 
pensator. We continue to assume Hypothesis 9.1. 


We will denote the resolvent set of [A-BF+F (C^f-C)] by 
p ( [A-BF+F (C q f-C ) ] ) . 

Theorem £.5. There exists a real number ^ suc h that, if Re(s) > a x , then s e 
p([A n -fi n F n +F n (CQF n -C n )] ) for all n, and ® n (s) converges to 5(s), uniformly in 
compact subsets of such s. 


Proof . The operator [A-BF+F (C^f-C)] is obtained from a contraction semigroup 
generator by perturbation with bounded operators, and the approximations to 
the perturbing operators are bounded in n, by strong convergence. In view of 
this, close examination of the basic approximation scheme in Section 4.1 will 
show that there exists a bound of the form M^exp(a^t), independent of n, for 
the semigroups generated by [A n ~B n F n +^ n (CQF n -C n )l . Also, these semigroups 
converge strongly to the semigroup generated by [A-BF+£(C q f-C)] , according to 
[G3, Theorem 6.6]. For Re(s) > a^ then, the resolvent operator in ® n (s) con- 
verges strongly to that in 5(s), uniformly in compact s-subsets, by [Kl, page 
504, Theorem 2.16, and page 427, Theorem 1.2]. □ 
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This result leaves much to be desired. For example, it does not guaran- 
tee that any subset of the imaginary axis will lie in 

p( [A -b F +F (C A F - CJ]) for sufficiently large n, even if all of the ima- 
nnnnunn 

A 

ginary axis lies in p( [A-BF+F(CqF-C)] ) • As with the convergence of the 
closed-loop systems, we can get more for certain important cases. 

Remark 9.6. If the open-loop semigroup T(*) (whose generator is A) is an ana- 
lytic semigroup, then there exist real numbers a, 9 and M, with 0 and M posi- 
tive, such that p( [A-BF+F(CqF-C)] ) contains the sector (s : larg(s-a) | < 2 + ^ ' 
and for each s in this sector, 

I I (si- [A-BF+F (C 0 F-C)] ) _1 II i M/|s-a| (9 17 ) 

Theorem 9.7. i) If the basis vectors of the approximation scheme are the 
natural modes of undamped free vibration and the structural damping does not 
couple the modes, then each s in p( [A-BF+F (C^f-C)] ) is in 

p([A n -B n F n +F n (C 0 F n _C n )] ) for n sufficiently large and « n (s> converges to S(s) 

A 

as n — > ® , uniformly in compact subsets of p( [A-BF+F(C qF-C) 1 ) . ii) If, 
additionally, T(») is an analytic semigroup, then ® n (s) converges to S(s) uni- 
formly in the sector described in Remark 9.6. 


Proof , i) In this case, we have also 


F n p En^ sI %iComp* F n* 


(9.18) 

where ^ n Q om p is the operator on D(A) defined by (9.11). The result follows 

from (9.8) and (9.9) and the fact that ^ nCoai p converges in norm to 

[A-BF+F(C 0 f-C) 1 . ii) The result follows from i) and a bound on 

(si- A „ ) -1 for large |s| that is obtained from the Neumann series in view 

nComp' 1 1 
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of (9.17) and the uniform-norm convergence of 

Theorem 9.8. If A has compact resolvent, then ^ n (s) converges to $(s) for 
each s e p( [A-BF+F(C qF-C)] ) , uniformly in compact subsets. 

Proof . As a result of Theorem 4.4, the resolvent operator in ® n (s) converges 
in norm to the resolvent operator in S(s) for sufficiently large real s. 

After an artificial extension of A^ to then, the present theorem follows 
from [Kl, pages 206-207, Theorem 2.25]. 
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10. Closing the Loop in the Example 


As in Definition 7.3, the optimal closed-loop system is formed with the 
optimal infinite dimensional compensator, which consists of the optimal con- 
trol law for the distributed model of the structure applied to the output of 
an optimal infinite dimensional state estimator. This optimal control law is 
the limit of the approximating finite dimensicnal control laws in Section 6. 

In this section, we first approximate the infinite dimensional estimator, as 
in Section 8. and then apply the approximating control laws in Section 6 to 
the approximating finite dimensional estimators to produce a sequence of fin- 
ite dimensional compensators that approximate the optimal compensator. 

10.1 The Estimator Problem 

We assume that the only measurement is the rigid-body angle 0 and that 

A _4 

this measurement has zero-mean Gaussian white noise with variance R 
We model the process noise as a zero-mean Gaussian white disturbance that has 
a component distributed uniformly over the beam, as well as two concentrated 
components that exert a force on the tip mass and a moment on the hub^ For 
this disturbance, the covariance operator Q has the form (7.23) with Q Q = 0, 

= 0 and $2 = 

We construct the approximating estimators as in Section 8.1. The gain 
for the n th estimator is given by (8.13) with the solution to the Riccati 
matrix equation (8.14). For the rigid-body measurement, the matrix C n is 
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According to (8.15), the matrix Q n is 



since W n is the matrix in (4.33). 
matrix.) Recall from Section 6.3 
elements. 


( 10 . 2 ) 

(As always, M~ n is the inverse of the mass 
that n = 2n g + i where ^ ls the number of 


Our only use for the functional estimator gains is to measure the conver- 
gence of the finite dimensional estimators to the optimal infinite dimensional 
estimator. To see the convergence of the approximating estimator gains, we 
compute the approximating functional estimator gains as in Section 8.3. Like 
the functional control gains, the functional estimator gains have the form 

f 0^,) , 

(10.3a) 

8 = (<l g'*gJg) ' 

(10.3b) 

and the corresponding approximations have the form 


A 

= , B ) . 

n fn'*gn' H gn' ' 

(10.4a) 

8 n = (<x gn’^gn’ p gn > • 

(10.4b) 

Mernxk 10 .1 We cannot guarantee as much about convergence for the approximat- 
ing estimators as we could for the approximating control problems in Section 
6. Since the damping in this example does not couple the natural modes and 
the rigid-body mode is observable, we would have Part i) of Theorem 8.5 if we 
were using the natural mode shapes as basis vectors. Therefore, we know at 
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least that a solution to the infinite dimensional Riccati equation (7.16) 
exists and that the infinite dimensional estimator that we want to approximate 
exists. The numerical results indicate that the solutions to the finite 
dimensional Riccati equations are bounded in n and that the functional esti 
mator gains converge in norm. The rigid-body mode prevents our guaranteeing a 
priori all the convergence that we want. If a torsional spring and damper 
were attached to the hub in the current example, we would have coercive stiff- 
ness and damping and Theorem 8.5 ii) would guarantee that the solutions to the 
finite dimensional Riccati equations converge strongly and that the functional 
estimator gains converge in norm for the basis vectors used here. Also, see 
Remark 6.2 and Remark 8.6. □ 

For damping coefficient c Q = 10 -4 , Figures 10.1 and 10.2 show 0 fn and and 

Tables 10.1 and 10.2 list the the scalars a fn ,o gn and P gn - Since 0 fn <O> 

(o) = 0, the convergence of implies the convergence of P fn = 

as in the control problem, P fn is not an independent piece of information 
about the estimator gains while, as far as our results gp. P gn is. We main- 
tain analogy with the control problem and list only P gn in the table. 
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'o.oo ib. oo 2b. 

Eisuce 10.1a. 
Damping 


) 30.00 <10.00 50.00 sb.00 7b .00 ab.OO Sib. 00 100.00 

f I 

Functional Estimator Gain Component 
0 = 10~ 4 ; estimator R = 10 -4 


1 0V0O T&ToO 2V.00 30. 00 vb 00 50. 00 tbT 00 70. 00 80. 00 90.00 

Figure lO.IJ*. Functional Estimator Gain Component d 

*gn 

Damping c Q = io~ 4 ; estimator R = 10 -4 


iio 


number of elements n 


2 ,3 ,4 ,5 .8 


number of element n 


2 . 3 . 4 . 5 . 8 









n e 

“fn 

°gn 

Pgn 

2 

5.0358 

12.680 

-1334.9 

3 

5.2514 

13.789 

-1455.4 

4 

5.3195 

14.149 

-1495.7 

5 

5.3478 

14.300 

-1512.2 

8 

5.3611 

14.371 

-1520.1 


Tatifi lO.i. Scalar Components of Functional Estimator Gains 

Damping coefficient c Q = io -4 ; estimator R = 10 -4 

number of elements n = 2, 3, 4, 5, 8 

6 


“fn a gn Pj 


4 5.3195 14.149 -1495.7 

6 5.3567 14.347 -1517.5 

8 5.3611 14.371 -1520.1 

10 5.3623 14.377 -1520.8 

latifi 10.2. Scalar Components of Functional Estimator Gains 

Damping coefficient c Q = 10 -4 ; estimator R = 10 -4 

number of elements n = 4, 6, 8, 10 

V 
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Figures 10.3a and 10.3b and Table 10.3 give the numerical results for the fin- 
ite dimensional estimators when the structural damping is zero. While Theorem 
8.4 says that the solutions to the finite dimensional Riccati equations for 
these estimators will not converge when the damping is zero and Q is coercive 

A 

on E, we have no result to predict the convergence for zero damping when Q is 

not coercive (even though is coercive on H) . From the numerical results 

A 

though, f does not appear to converge. 


°fn 

a_„ 

gn 

Pgn 

5.0730 

12.868 

-13 54.4 

5.3390 

14.253 

-1506.0 

5.4417 

14.806 

-1568.0 

5.4894 

15 .067 

-15 96.3 

5 .5398 

15.345 

-1627 .2 


Table 10..3,. Scalar Components of Functional Estimator Gains 

_4 

Zero damping; estimator R = 10 
number n g = 2 , 3 , 4 , 5 , 8 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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10.2 Approximation of the Optimal Compensator 


Finally, for the damping c 0 = lO -4 , R = .05 in the control problem and R 
= 10 -4 in the estimator problem, we construct the finite dimensional compensa- 
tor in Figure 9.1; i.e, for each n = 2n g + 1. we apply the n£ h control law 
represented by the functional gains in Figure 6.4 and Table 6.4 to the output 

of the n fc k estimator represented by the functional gains in Figure 10.2 and 
e 

Table 10.2. As the number of elements increases, the transfer function m 
(9.16) of the finite dimensional compensator converges to the transfer func- 
tion in (7.14) of the optimal infinite dimensional compensator, as described 
in Section 9.3. Theorem 9.5 and Remark 9.6 apply. Figure 10.4 shows the fre- 
quency response (bode plots) of the finite dimensional compensators f or 4 , 6 . 

8 and 10 elements. The phase plot is for 10 elements only. These plots indi- 
cate that the finite dimensional compensator for eight or more elements is 
virtually identical to the optimal infinite dimensional compensator, as 
predicted by the functional gain convergence in Figures 6.4 and 10.2. 
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10.3 Comments on the Structure and Dimension of the 


Implementable Compensators 

Though this paper does not address the problem of obtaining the lowest- 
order compensator that closely approximates the infinite dimensional compensa- 
tor, we should note that the compensators based on eight and ten elements here 
are unnecessarily large because the finite element scheme that we chose is not 
nearly the most efficient in terms of the dimension required for convergence. 
(The dimension of the first-order differential equation in the compensator is 
2(2ne+l) .) We used cubic Hermite splines here to demonstrate that the finite 
element scheme most often used to approximate beams in other engineering 
applications can be used in approximating the optimal compensator. In [G5] . 
we compare the present scheme with one using cubic B-splines and one using the 
natural mode shapes as basis vectors. The natural mode shapes yield the 
fastest converging compensators, but the B-splines are almost as good. The 
only advantages of the Hermite splines result from the fact that the coding to 
build the basic matrices (mass, stiffness, etc.) is simpler than for B splines 
and the fact that, before the Riccati equations based on, say, ten natural 
modes are solved, a much larger finite element approximation of the structure 
must be used to get the ten modes accurately. 

To understand the redundancy in the large finite dimensional compensators 
here, it helps to consider the structure of the optimal compensator. It is 
based on an infinite dimensional state estimator that has a representation of 
each of the structure's modes. In the present example, the optimal compensa- 
tor estimates and controls the the first six modes significantly, the next 
three modes slightly, and virtually ignores the rest. This observation is 
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based on the projections of the functional gains onto the natural modes and on 
comparison of the open- loop and closed-loop eigenvalues. (See [G5] for more 

detail, including the spectrum of closed-loop system — which is stable 

obtained with the ten-element compensator here.) The infinite dimensional 
compensator then has an infinite number of modes that contribute nothing to 
the input-output map of the compensator. These inactive modes are just copies 
of all the open-loop modes past the first nine. They can be truncated from 
the compensator without affecting the closed-loop system response signifi- 
cantly. The number of active modes in the compensator — i.e. , the modes that 
contribute to the input-output map — depends on the structural damping and 
the Q's and R's in the LQG problem statement. (See the discussion in Section 
6.4 about the effect of damping and control weighting on performance.) 

The compensator computed here based on ten elements has 21 modes 
(although we did not do the computations in modal coordinates). Nine of these 
compensator modes are virtually identical to the nine active modes in the 
infinite dimensional compensator, and the twelve inactive modes are approxima- 
tions to the tenth through twenty-first open-loop modes of the structure. The 
inactive modes result from the large number of elements needed to approximate 
the active compensator modes accurately. Now that we essentially have the 
optimal compensator in the ten- element compensator, we could truncate the 
twelve inactive modes and implement a compensator with nine modes. And we 
probably could reduce the compensator even further using an order reduction 
method like balanced realizations. 
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11. Conclusions 


For the deterministic linear-quadratic optimal regulator problem for a 
flexible structure with bounded input operator (the B Q i n (2.1)), the approxi- 
mation theory in Sections 4 and 5 is reasonably complete. The most important 
extensions should be to the corresponding (very difficult) problem with 
unbounded input operator, for which there exists little approximation theory. 
Because of the different kinds of boundary input operators, stiffness opera- 
tors and structural damping, all of which must be considered in detail when B Q 
is unbounded, it seems unlikely that the approximation theory for the 
unbounded-input case can be made as complete as the theory here. 

The convergence results in Section 8 for the estimation problem are less 
complete than those for the control problem because rigid-body modes present 
more technical difficulties for the proofs in the estimator case. However, 
our analysis and numerical experience suggest that the difficulties only make 
the proofs harder and that the convergence in the estimation problem is ident- 
ical to the convergence in the control problem, and that controllable and 
observable rigid-body modes make no qualitative difference in either problem. 

Where we would most like substantial improvement over the results of this 
paper is in Section 9.2, which considers how the approximating closed-loop 
systems obtained by controlling the distributed model of the structure with 
the finite dimensional compensators converge to the optimal closed-loop sys- 
tem, obtained with the infinite dimensional compensator. Theorem 9.4 gives us 
what we want for problems where the damping does not couple the natural modes 
of free vibration and the natural mode shapes are the basis vectors for the 
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approximation scheme. In particular, this theorem says that, if the optimal 
closed-loop system is uniformly exponentially stable, then so are the approxi- 
mating closed-loop systems for sufficiently large order of approximation. We 
have verified numerically the stability of the approximating close-loop sys- 
terns for the example in Sections 6 and 10, where the basis vectors are not the 
modes. This example and others have made us suspect that Theorem 9.4 is true 
when the basis vectors satisfy Hypothesis 4.1 only and when the damping cou- 
pies the modes. 

Another possible approach to analyzing the convergence of the approximat- 
ing closed-loop systems to the optimal closed-loop system is to use the 
input-output description in frequency domain. Results like those in Section 
9.3 are useful for this, although for the closed-loop stability we want, we 
probably need the transfer functions of the finite dimensional compensators to 
converge more uniformly on the compensator resolvent set than we have proved 
here. In our example. Figure 10.4 indicates that these transfer functions 
converge uniformly on the imaginary axis, but we have no theorem that guaran- 
tees this. 
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APPENDIX 


Errata for [Gl] 

In the first paragraph of the proof of Theorem 2.1 on page 
689 of [Gl], the first sentence should be: 

If a dissipative operator is invertible, its inverse is 
dissipative . 

At the beginning of the fifth line of the same paragraph, the 
expression (otx+y) should be deleted the first time it occurs. 

The next-to-last sentence of the paragraph should be: 

Hence, if a densely defined maximal dissipative operator has 
dense range, its inverse is maximal dissipative. 

The theorem is correct as stated. 


In the current paper, we use Theorem 2.1 of [Gl] to conclude 
that the operator A defined in Section 2 is maximal dissipative (see 

o 

(2 . 10) - (2 . 12 )) and that the operator A in (2.16) has a unique maximal 
dissipative extension. 
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